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Abstract 


This  paper  presents  an  improvement  for 
the  mean-stream  line  method  -  using  the 
analytical  solution.  This  method  can  drasti¬ 
cally  speed  up  the  calculation  and  increase 
the  accuracy  of  the  calculation.  In  this 
paper  it  presents  simple  equations  and  tables 
of  functions  for  convenience  of  calculation 
as  well  as  examples  of  such  calculation.  In 
order  to  convenience  the  designing,  this  paper 
also  analyzes  preliminarily  as  to  how  the 
design  parameters  should  be  chosen  in  order 
to  obtain  the  desired  turbine  blade  shape. 

Using  these  methods  and  related  information, 
it  can  be  done  in  one  day  to  complete  the 
design  of  a  good  turbine  cascade  circulated 
by  compressible  flow  (if  only  calculation  of 
Incompressible  flow  is  necessary,  the  speed 
of  computation  can  still  be  rised  several 
time).  This  paper  finally  also  gives  a 
comparison  between  its  results  and  the 
experimental  data  and  presents  several  possible 
directions  under  which  the  mean-stream  line 
method  can  still  be  further  developed. 


SYMBOLS 


A  -  trigonometric  function  of  6 

a  -  width  of  the  exit  in  the  throat  region  in  the  duct 
between  blades 

*This  paper  has  been  presented  in  the  First  National  Meeting 
of  Engineering  Thermal  Physics  in  August,  1965. 
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trigonometric  function  of  6 ,  or  the  width  of  the 
blade  in  the  Z  direction 

arc  length  of  the  blade 

trigonometric  function  of  g 
isobaric  specific  heat  of  the  gas 

diameter  of  the  inner  circle  of  the  duct  between  blades 
trigonometric  function  of  3 
trigonometric  function  of  g 
enthalpy  of  the  gas 

blade  surface  or  the  arc  length  of  the  blade  surface 
Mach  number  of  the  gas  flow 

the  distance  in  z-direetion  of  the  front  fringe  of  the 
blade  at  the  point  where  the  thickness  in  the  y-direction 
is  the  maximum 

trigonometric  function  of  g 

pressure  of  the  gas 

an  arbitrary  physical  quantity 

gas  constant 

radius 

trigonometric  function  of  g 
cascade  distance  of  blade  cascade 
gas  flow  velocity 

coordinate  direction  (direction  of  the  forehead  line 
of  the  cascade) 

coordinate  direction  (axial  direction  of  the  turbine 
machinery ) 

the  absolute  value  of  the  distance  in  the  y  direction 
between  the  mean-stream  line  and  the  surface  of  the 
blade 

an  angle  (starting  from  the  z-coordinate  and  giving 
counter  clockwide) 

the  thickness  of  the  blade  in  the  y-direction  at  any 
point 

ratio  of  specific  heats  of  gas 
a  constant 

a  constant;  or  the  gas  flow  velocity  divided  by  the 
critical  velocity,  i.e.  velocity  coefficient 

a  constant 
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p  -  the  density  of  gas 
w  -  a  function 

SUPERSCRIPTS 

*  -  relative  value  (divided  by  the  inlet  parameter) 

'  -  derivative  with  respect  to  z 

.  -  parameter  at  the  rear  forehead  line  of  the  cascade 

SUBSCRIPTS 

o  -  stagnation  value 

1,2, 

3. . .n  -  different  functions  or  constants 
B  -  incompressible 

e  -  far  away  downstream  of  the  cascade 

i  -  far  away  upstream  of  the  cascade 

m  -  values  on  the  mean-stream  line  or  values  on  any 
selected  stream  line 

max  -  maximum  value 

p  -  Inner  arc  of  the  blade 

s  -  back  arc  of  the  blade 

y  -  y-direction  component 

z  -  z-direction  component 

60?5  -  values  on  the  stream  line  which  has  a  60%  flow  of  that 
at  the  back  arc. 

I.  Mean-Stream  Line  Method 

The  condition  of  the  gas  flow  circulation  through  the 
planar  cascade  has  an  extremely  large  effect  on  the  character¬ 
istics  of  the  turbine  machinery.  Therefore  the  calculation 
of  the  condition  of  gas  flow  circulation  for  a  given  planar 
cascade  (Forward  Problem)  and  the  designing  of  a  planar  cas¬ 
cade  which  has  good  aerodynamic  properties  and  satisfies  cer¬ 
tain  actual  requirements  (reverse  problem)  are  both  important 


subjects  to  study  for  people  working  in  the  turbine  machinery. 

Up  to  date,  there  has  been  a  lot  of  work  done  in  this  area. 

There  are  various  methods  available.  Among  them  there  are 
methods  which  are  based  on  relatively  rigorous  theoretical 
basis  such  as  the  conformal  transformation  method  and  so  on. 
However,  the  actual  calculation  load  is  relatively  large. 

In  addition,  it  is  very  difficult  to  extend  to  compressible 
flow  and  flow  on  an  arbitrarily  revolving  surface.  Pure 
numerical  solution  has  been  used  to  consider  the  compressibility 
of  gases  but  the  computational  load  is  too  large  and  it 
casually  requires  an  electronic  computer  to  carry  out  the  com¬ 
putation.  Therefore,  on  the  other  '.jnj,  some  relatively  simple 
approximation  methods  energed .  For  example,  the  mean-stream 
line  method^1^  proposed  by  Professor  Wu  Ching-hua  is  one  of 
the  effective  methods  to  calculate  planar  cascades.  It  has  a 
relatively  sophisticated  theoretical  basis  to  easily  consider 
the  compressibilities  of  gases.  Its  computational  load  is  not 
large  and  results  with  sufficient  accuracy  can  be  obtained 
rapidly.  In  addition,  the  thickness  distribution  required  in 
the  design  work  can  be  easily  found  in  order  to  assure  the 
strength,  rigidity,  and  cooling  requirements.  Besides,  it 
can  be  extended  to  the  more  generalized  conditions  of  flow  on 
an  arbitrary  revolving  surface.  Therefore,  it  has  already 
become  a  generally  used  computational  method  worldwide. 

The  essence  of  the  mean  stream  line  method  is  as  follows: 

If  the  shape  of  a  stream  line  inside  the  duct  of  the  cascade  is 

known  (such  as  the  mean  stream  line  y  =  y(z)  which  equally 

m 

divides  the  flow  on  the  duct)  and  the  variation  pattern  of  any 
gas  flow  parameter  on  it  [such  as  («(',),*«(*)  ]  (refer  to 

Figure  1),  then  it  is  possible  to  use  the  steady,  isoentropic, 
continuous,  non-rotational,  and  ideal  gas  conditions  in  aero¬ 
dynamics  to  derive  the  partial  derivatives  of  various  orders 
of  the  gas  parameters  on  that  stream  line  with  respect  to  the 
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y  direction.  The  actually  applicable  equations  are  as 
follows^1^: 
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The  q  in  equation  (9)  can  represent  any  one  of  the  gas 
flow  parameters.  In  the  actual  calculations  of  the  turbine 
cascade,  only  the  first  three  terms  of  the  Taylor  series 
expansion  will  provide  sufficient  accuracy. 

In  order  to  solve  the  ’’Forward  Problem",  at  the  beginning 
of  the  calculation  it  is  necessary  to  estimate  ym  =  y(z) 
and  (pik,)„  =o(0.  under  the  given  duct  condition. 
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In  Reference  [1],  it  was  believed  that  the  center  line 
of  the  duct  and  be  used  as  the  first  approximation  of  the 
mean-stream  line  (the  amount  of  flow  equals  on  both  sides  of 
this  line)  and  the  variation  of  the  parameter  PII'«  along  the 
mean  stream  line  is  inversely  proportional  to  the  y-direction 
width  of  the  duct  (i  .e .  (pn\) On  this  basis  of 
these  a  preliminary  calculation  can  be  carried  out.  As  for 
whether  its  results  are  correct,  we  can  integrate  from  z’s  on 
the  mean  stream  line  along  the  y  direction  on  both  sides  to 
determine  whether  the  obtained  flow  agrees  with  the  incoming 
flow.  Reference  [1]  also  gave  a  modification  method  when 
disagreement  was  found.  Later,  in  Reference  [3],  it  presented 
a  method  to  more  accurately  estimate  the  position  of  the  mean 
streamline . 

For  the  "Reverse  Problem",  if  >«.=  ?(*)  and  (PH',)m=  o>  (  * ) 
are  given,  then  all  the  flow  parameters  in  the  flow  field  can 
of  course  be  obtained.  In  the  meantime,  the  boundary  of  the 
duct  (blade)  can  also  be  defined  based  on  the  concept  of 
equilibrium  flow.  Naturally,  to  obtain  good  cascades,  .’-,=  >’(«) 
and  (pu\)„=  a ( z  )  also  must  have  certain  limitations  and  matching 
relations.  Reference  [3]  introduced  a  method  using  past  exper¬ 
ience  to  plot  the  shape  of  the  blade  and  to  obtain  the 
approximate  variation  regularity  between  the  mean  stream  line 
and  the  duct  width.  Afterwards,  the  preliminary  calculation 
is  carried  out  followed  by  gradual  improvements. 

II.  One  Way  to  Increase  the  Computational 
Speed  and  Accuracy 

In  the  method  of  using  existing  mean  stream  line  to 
design  cascades  (the  Reverse  Problem),  as  described  above, 
because  usually  the  mean  stream  line  and  the  data  of  one  para¬ 
meter  at  several  points  on  it  are  given  (in  the  past 


was  obtained  from  consideration  of  the-  Made  thickness), 
therefore  the  numerical  differentiation  method  is  used  to 
calculate  the  derivatives  ~f.~  along  the  stream  line  in 
equations  (!i)  ~  (7d).  This  part  of  the  calculation  is  more 
time  consuming  and  is  difficult  to  obtain  very  nigh  accuracy. 

It  is  also  easier  to  make  mistakes.  Because  of  the  use  of 
numerical  differentiation  equations,  the  error  of  one  mis¬ 
takenly  calculated  point  will  affect  the  accuracy  of  the 
numerical  values  thereafter.  References  [1]  ~  [3]  suggested 
the  plotting  of  curves  in  a  table  after  these  derivatives  are 
calculated  in  order  to  discover  mistakes  early  and  to  "polish" 
the  data. 

The  mean  stream  line  method  does  not  require  many  com¬ 
putational  terms  which  is  one  of  its  advantages.  But  at  least 
It  times  of  numerical  differentiation  must  be  calculated.  If 
the  five  point  differentiation  eauation^-1^  is  used,  then  the 
computation  of  each  point  involves  5  multiplications,  1  addi¬ 
tion,  and  1  division.  In  addition,  in  order  to  obtain  a 
certain  significant  figures  with  accuracy,  it  is  necessary  to 
use  more  significant  figures  in  the  beginning  of  the  calcula¬ 
tion.  This  then  becomes  a  barrier  which  may  not  necessarily 
exist  against  increasing  the  computational  speed.  Besides,  the 
polishing  in  the  numerical  calculation,  although  is  an  effec¬ 
tive  method,  yet  it  does  not  absolutely  assure  the  accuracy 
of  the  calculation.  Therefore,  it  is  natural  to  have  the 
following  thoughts:  Is  it  possible  to  change  the  design  of 
the  mean  stream  line  method  from  variation  of  parameters  to 
the  analytical  form?  This  allows  the  use  of  the  analytical 
method  in  the  differentiations  and  the  rearrangement  afterwards 
so  that  the  complicated  calculations  described  above  can  be 
eliminated.  By  doing  so  it  is  not  necessary  to  gradually  carry 
out  calculations  of  numerical  differentiation  to  reduce  the 
possibility  of  computational  error.  It  is  also  not  necessary 
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to  publish  the  computational  results  and  can  obtain  the  deri¬ 
vatives  to  an  arbitrary  degree  of  accuracy  desired.  For 
the  logical  point  of  view,  this  method  has  no  difficulty  in 
the  "Reverse  Problem".  As  long  as  we  can  find  the  proper  mean 
stream  line  and  the  functional  relationship  between  the  flow 
parameter  on  it  and  the  z-coordinate .  These  are  not  difficult 
to  obtain  from  the  existing  cascade  data  and  the  computational 
experience  acquired  later.  Of  course,  if  the  equations 
expressed  by  the  analytical  form  are  too  complicated  then  its 
meaning  is  no  longer  significant.  But  from  the  following 
section,  it  can  be  seen  that,  to  the  contrary,  the  use  of 
analytical  relations  can  simplify  the  equation  itself  which  is 
more  advantageous  in  computation.  Let  us  elaborate  this  analy¬ 
tical  method  in  the  following. 

III.  Analytical  Solution  of  the  Reverse  Problem 


1.  Incompressible  Flow 


Let  us  first  take  a  look  at  the  special  solution  to  the 
incompressible  flow.  At  this  time  equations  (1),  (4),  (5), 
(7),  and  (7a)  become  the  following,  respectively.  (The  para¬ 
meters  in  this  section  and  in  Table  2  in  the  Appendix  are 
with  respect  to  the  mean  stream  line.  For  convenience,  the 
subscript  is  omitted  with  exception  of  those  specified  in 
the  equations): 


«-(+).-• 

(10) 

«H-(T  'I- ).“»’» 

(ID 

(12) 
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If  we  choose  the  analytical  equation  of  the  mean  stream 
line  to  be  y  =  v(z)  .  In  addition,  we  also  choose  W=W(z) 
as  the  varying  parameter  of  design  instead  of  M'.swCz). 

This  is  apparently  possible  in  solving  the  reverse  problem  and 
offers  convenience  in  computation.  It  makes  easier  to  control 
the  surface  velocity  distribution  of  the  obtained  blade  shape. 
Actually  after  =  and  y  =  v  ( z )  are  given,  then  in 

principle  it  is  completely  possible  to  transform  into  VJ-W(z). 
Thus  using  equations  (2),  (3),  and  the  above  four  equations 
as  well  as  tan  P  =  y ' f 3  )  ,  it  is  possible  to  derive  the  equation 
of  the  analytical  solution.  However,  the  operation  is  rather 
complicated  and  the  resulted  equation  is  more  difficult  to 
use.  Therefore,  it  is  better  to  obtain  8  “taiT'O'C  z  )  )  =  8  (  *  ) 
from  y  =  y(z)  as  the  varying  parameter  and  then  differentiate 
W  and  W  with  respect  to  z  to  carry  out  the  derivation  opera- 

y  ^ 

tion.  It  is  then  possible  to  get  the  simple  equations  of  the 

partial  derivatives  of  W  and  W  in  the  y  direction. 

y  z 


(&W  \ 

(15) 

+  A.'  +  A? ■  W ’ 8'  ~  Bt •  W  (1 6) 

.  (rar),  (17) 

+  + Bf\V '  (18) 


where  A^,  Bg  are  the  simple  trigonometric  functions  of 

8.  For  example: 


^  =  0064*8  sin  p  (12  cos "8  —  5  ) 


(19) 
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/f,=  -  cos  P  (16  cos‘P-  14  cosJp  ^  1  ) 

(20) 

B,  =  cos  P  (  2  coerp  -  1  ) 

(21) 

Bj=  2  cos:p  sin  P 

(22) 

The  remaining  equations  of  B  can  be  found  in  Table  3  in  the 
Appendix . 

But  in  the  actual  application,  it  is  more  convenient  to 
calculate  after  the  values  of  the  partial  derivatives  of  W- 
and  W  are  obtained.  It  is  because  that  the  former  can  be  used 
to  determine  the  boundaries  of  the  duct  and  the  latter  can  be 
used  to  directly  obtain  the  velocity  distribution  with  com¬ 
bining  the  two  velocity  components  to  get  the  velocity.  The 
partial  derivative  value  of  W  can  be  obtained  as  follows: 

Since 

i  r; 


Let  us  carry  out  first  order  and  second  order  partial  differ¬ 
entiations  on  both  sides  of  the  above  equation.  After 
rearrangement  and  simplification  we  get 
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then  substituting  equations  (15)  ~  (18)  into  the  two  above 
equations,  we  get 


(*£).“»»•»'< V+BJV  (25) 


=  B,-!rP'5  +  2Bi0.|{’P'  +  BJi-ir'P'  +B;,ir'  +  b„- 
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where  ^  are  all  simple  trigonometric  functions  of  B. 


The  actual  equations  can  be  found  in  Table  3  in  the  Appendix. 


From  these  it  is  found  that  for  incompressible  flow  the 
first  order  and  second  order  partial  derivatives  in  the  y 
direction  on  the  mean  stream  line  are  simple  linear  function 
of  the  combination  of  products  of  W,  6  and  their  partial 
derivatives  with  respect  to  z.  Because  the  coefficients  in 
front  of  each  term  are  all  known  simple  trigonometric  functions 
of  $,  table  (see  Table  3  in  the  Appendix)  can  be  prepared  in 
advance  for  future  use.  After  doing  so  it  is  very  quick  to 
calculate  the  partial  derivatives  in  the  y-directions  on  the 
mean  stream  line.  For  incompressible  flow,  after  the  familiar¬ 
ization  of  the  process,  it  takes  about  or  less  than  two  hours 
to  calculate  over  ten  points  of  a  duct. 


In  the  next  step  of  determining  the  position  of  the 
boundary  of  the  duct,  it  is  possible  to  use  the  partial  deri¬ 
vative  values  everywhere  on  the  mean  stream  line  and  the  con¬ 
tinuity  equation  to  derive  a  third  order  algebraic  equation  to 
replace  the  numerical  integration  method  used  in  the  past. 


Since 


0 


and 


can  get 


2  ( V )  Ay> + = -r»V 


Similarly 


-  -r»v 


*The  W  Tn  these  two  equations  represents  the  axial  velocity 
component  in  the  direction  of  gas  flow  at  any  point  in  the 
duct.  It  is  different  from  the  axial  velocity  component  on  the 
mean  stream  line  in  other  places  of  this  section. 
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where  absolute  values  are  used  for  Ay,  and 


Of  course,  it  is  not  a  simple  matter  to  rigorously  solve 

the  third  order  algebraic  equation  (although  it  can  be 

*  * 

obtained  completely).  Eut  its  approximate  solution  with 

a  certain  accuracy  can  be  easily  obtained  graphically  or  using 
the  interpolation  method.  Based  on  the  experience  in  calcula¬ 
tion,  for  blade  of  ordinary  size  (arc  length  in  the  order  of 
several  tens  mm),  it  is  only  necessary  to  suostitute  the  two 
integer  value  in  millimeters  of  Ay  near  the  solution  (the 
approximate  value  can  be  estimated  easily  from  the  possible 
shape  of  the  blade)  into  equation  (27)  or  (28)  and  then  use 
the  linear  interpolation  method  to  obtain  the  solution  which 
is  pretty  accurate.  After  the  boundary  of  the  duct  is  known, 

we  can  use  equation  (9)  and  W  as  well  as  its  partial  deriva- 

m 

tive  value  in  the  y-direction  to  obtain  the  velocity  distri¬ 
bution  on  the  surface  of  the  blade. 

2.  Compressible  Plow 


For  compressible  flow,  it  is  also  completely  possible 
to  use  a  similar  method  as  described  above  to  obtain  the 
corresponding  equations.  As  discussed  before,  at  this  time 
we  should  find  the  partial  derivatives  of  Pit’,  and  W  in  the 
y-direction  and  not  those  of  and  Wz  in  order  to  make 
further  calculation  more  convenient.  Equations  (23)  and  (24) 
are  still  valid  at  this  time.  We  also  noticed  that: 
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**  Recently,  in  Reference  [7],  it  presented  a  method  using 
the  concept  of  using  the  flow  function  to  directly  calculate 
Ayp  and  Ayg. 

-a-?" 


.  ** 


Then  from  equations  (1)  ~(8)  after  complicated  operation,  we 
can  obtain  the  following  results: 
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where  every  B  is  also  a  trigonometric  function  of  B.  The 
actual  equations  can  be  found  in  Table  3  in  the  Appendix. 


From  here  it  can  be  found  that,  for  compressible  flow 
although  there  are  so  many  terms  yet  it  is  also  the  addition 
of  some  simple  products.  The  coefficients  of  each  term  can  be 
tabulated  in  advance  (see  Table  3  in  the  Appendix)  for  use 
later.  Therefore,  the  computation  is  also  relatively  quick 
and  it  takes  a  day  or  so  to  complete.  The  determination  of 
the  position  of  the  boundary  of  the  duct  at  this  time  can  also 
use  the  third  order  algebraic  equation  and  not  the  numerical 
integral  method: 
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When  y  =  y(»)  is  given,  the  values  of  6,  6’  and  6"  can 
be  obtained  from  the  following  equations: 


P  *=tan_,Cy'  (  *  )) 
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(36) 

(37) 


P'  =  ~  i  *  y>3-  =  =  B,y* 

p'  =  cos:3(y'-y':.2  ccw  P  sin  p  )  =  B:y'  -  2/?lby’1 


Therefore,  when  some  ordinary  geometric  curves  are  used  as  the 
mean  stream  line,  the  equations  of  each  y-direction  derivative 
can  also  be  simplified  using  the  above  three  equations.  Table 
2  in  the  Appendix  lists  the  equations  under  the  special  con¬ 
ditions  that  the  mean  stream  lines  are  six  different  curves, 
respectively.  The  c,-c„  Tj-r*  and  Hj-h,  etc.  in  the 

equations  ae  also  known  trigonometric  functions  of  6.  The 
expressions  and  values  can  be  looked  up  in  Table  3. 

IV.  Application  Examples 

In  this  section,  we  are  using  a  high  subsonic  pulsed  cas¬ 
cade  as  an  example  whose  design  requires  relatively  compli¬ 
cated  computations.  There  are  many  terms  to  be  omitted  to 
simplify  the  calculation. 

The  original  data  of  the  design  are  as  follows: 

Gas  Paramaters  before  the  Cascade: 

stagnation  temperature  r,,  =  324‘K  , 

stagnation  density  Po<  =  l.565  kg/m^, 
velocity  w'i-iso  m/sec 
angle  F.  =  5i.2*  ; 

Gas  Parameters  behind  the  Cascade: 

velocity  H'#  =  264  m/sec 
angle  P«=-60’  j 

The  working  medium  is  an  ideal  gas  with  Y=I.4 


t 

J 


The  geometric  shape  of  the  mean  stream  line  selected  in 
the  design:  when  P„>0  ,  it  is  a  circle  with  radius  r  =  50 

mm.  When  P,n<0  ,  it  is  a  suspended  chain  wave  v  =  'l cosh  (Q-)  +  X 

(mm)  where  *1  =  -30.77  mm,  0  =0.0255  1/mm,  and  X  =80.7;  mm 

When  the  combination  of  the  two  types  of  geometric  curves 
is  used  as  the  mean  stream  line,  it  is  necessary  to  make  sure 
that  the  values  of  and  are  the  same  on  both 

sides  of  the  junction  of  the  two  types  of  curves.  Only  by 
doing  so,  the  results  obtained  from  the  above  equations  will 
be  continuous  (when  only  the  first  three  terms  of  the  Taylor 
series  are  used).  For  ordinary  simple  geometric  curves,  this 
requirement  can  not  be  met  at  any  point. 

The  front  and  back  forehead  lines  of  the  cascade  are 
selected  right  at  the  places  where  the  Bm  of  the  mean  stream 
line  is  equal  to  the  given  and  Bg .  For  preliminary  design, 
this  selection  is  the  only  natural  choice.  But  from  the  dis¬ 
cussion  in  the  following  section,  we  can  see  that  this  choice 
should  be  recommended  to  be  used  in  general. 

The  velocity  distribution  of  the  mean  stream  line  in  this 
example  is  selected  as  follows:  When  Pm>0  J  use  W  =  constant 
=  200  m/sec.  When  Pm<0  ,  then  let  the  velocity  increase 
following  a  cosine  curve  to  the  exit  value  (270  m/sec),  i.e. 

! B'  =  235  —  35  cos^-~  re )  where  zg  is  the  value  of  z  coordinate  of 
the  back  foreward  line  of  the  cascade.  Although  at  P = 0  W 
and  W’  are  continuous,  W"  is  not  continuous.  The  continuities 
of  W,  W' ,  and  W"  all  should  be  required  at  the  junction  point 
of  various  velocity  variations.  But  it  is  generally  not  easy 
to  achieve.  Therefore,  for  the  velocity  distribution  on  the 
mean  stream  line  in  this  example,  we  can  imagine  the  following 
When  0  ,  use  a  constant.  When  ,  use  the  cosine 

curve  as  discussed  above.  As  for  the  segment  in  the  middle, 
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Table  1  (continued) 


Key  to  Table  1:  2.  Item,  3.  z(mm),  4.  y(mm),  5.  w(m/sec), 

C  .  .mm ) ,  7.  W"(M/sce  mm2 ) ,  8.  V'  (n/sec),  9-  equation, 

z 

10.  sc-e  equation  (3),  11.  equation,  12.  Based  on  6th  column 
of  t he  3rd  line  in  Table  2,  13.  Based  on  6th  column  of  the 
2nd  line  in  Table  2,  l1*.  see  equation  (3),  15.  item,  16.  s  (mm) 
17.  y  (mm),  18.  V.’(m/sec ) ,  19.  V/'  (m/sec.mm),  20.  V.'"  (m/sec  .nun2 )  , 
21.  W  _  (m/see),  22.  p  (kc/m3),  23.  (km/m2  sec),  24.  (m/sec.mm), 
25.  (ffl/sec.mm),  26.  look  up  the  table  of  aerodynamic  functions 
based  on  A,  27.  look  up  the  table  of  aerodynamic  functions 
based  on  A,  23.  see  the  2nd  column  in  line  6  in  Table  2,  29. 
see  the  2nd  column  in  line  10  in  Table  2,  30.  look  up  the 
table  of  aerodynamic  functions  based  on  A,  31.  look  up  the 

tab  e  of  aerodynamic  functions  based  on  A,  32.  see  the  6th 

column  in  line  6  in  Table  2,  33.  see  the  6th  column  in  line 
10  in  Table ^ 2. 

34.  p(Kg/m^),  35.  pWz  (Kg/m  sec),  36.  (M/sec.mm),  37.  M/sec.mm 
38.  see  col.  2  of  line  7  in  Table  2,  39.  see  col.  2  of  line  11 
in  Table  2,  40.  see  col.  2  of  line  8  in  Table  2,  41.  see  col. 

2  of  line  9  in  Table  2,  42.  see  col.  2  of  line  12  in  Table  2, 

43.  see  col.  6  of  line  7  in  Table  2,  44. see  col.  6  of  line  11 

in  Table  2,  45.  see  col.  6  of  line  8  in  Table  2,  46.  see  coi. 

6  of  line  9  in  Table  2,  47.  see  col.  6  of  line  12  in  Table  2 
48.  use  the  first  3  terms  of  equation  (9),  49.  use  the  first 

3  terms  of  equation  (9),  50.  use  the  first  3  terms  of  equation 

(9),  51.  use  the  first  3  terms  of  equation  (9) 


a  higher  order  of  smooth  curve  is  used  as  the  transition. 

Because  the  calculation  points  are  always  finite  and  discrete, 
therefore  transition  condition  in  between  calculation 

points  can  have  some  decree  of  arbitrariness.*  By  how  much 
the  velocities,  at  the  entrance  and  exit  on  the  mean  stream 
line  are  greater  than  the  given  velocity  values  in  the  front 
and  the  tael:  of  the  cascade,  by  how  much  the  value  of  *  at 
the  corresponding  point  is  greater  than  the  given  values  in 
front  and  at  the  back  of  the  cascade,  also  determines  the 
extent  of  dullness  at  the  front  fringe  of  the  blade.  The 

M  number  in  this  example  is  large,  therefore,  a  sharper 
design  is  used  in  order  to  avoid  a  partial  supersonic  effect 
at  the  entrance. 

The  cascade  distance  in  this  example  is  56.4  mm  corre¬ 
sponding  to  a  cascade  density  of  1.8. 

The  computational  results  and  procedures  are  shown  in 
Table  1  and  Figures  2  and  3.  From  the  calculated  results,  it 
can  be  found  that  the  coordinates  of  points  on  the  surface 
of  the  blade  shape  obtained  using  the  analytical  solution  of 
the  mean  stream  line  method  are  completely  smooth.  It  is  much 
smoother  than  the  solution  obtained  using  the  numerical  method 
(for  example)  refer  to  Reference  [3]).  This  also  indicates 
that  the  accuracy  of  the  solution  is  correspondingly  higher. 
Naturally,  as  described  above,  the  speed  of  calculation  is  also 
much  faster. 

For  the  design  of  a  cascade,  the  blade  profile  obtained 
in  this  example  can  be  considered  to  be  satisfactory:  even 
when  the  M  number  near  the  exit  reaches  0.8,  there  is  no 
supersonic  region  in  the  entire  duct.  There  is  even  some  mar¬ 
gin  (the  maximum  M  number  on  the  back  arc  of  the  blade  profile 

*  The  discontinuity  of  w’  at  the  junction  still  has  some  effect 
on  the  smoothness  of  the  calculated  results.  For  example,  the 
variation  of  t  in  the  region  where  is  not  most  ideal. 

YP 

But  the  effect  is  small  and  the  surface  of  the  blade  is  still 
smooth. 

zr 


is  about  0.9*0-  The  velocity  distribution  on  the  entire  blade 
is  also  satisfactory.  The  variation  is  uniform  O1', ir.  "'-'T  . 
The  compression  gradient  in  the  compression  region  which  exists 
in  the  exit  section  of  the  back  arc  of  an  ordinary  cascade 
is  not  large.  This  indicates  that  the  use  of  the  method  pre¬ 
sented  in  this  paper  to  design  good  planar  cascades  is  com¬ 
pletely  possible. 


•31 3 

Figure  2.  The  schematic  of  the  blade  profile  obtained  from 
the  sample  calculation  B1  -  51.2°,  B0  =  -60°,  =  0.8. 

Key:  1.  mean  line  of  the  inner  contact  circle  inside  the 
duct,  2.  mean  stream  line  ym  =  y(z),  3-  coordination  points 

obtained  from  calculation. 

(*•  =  !.«.  5  3- =0.2^.  a=0.30«) 

As  for  the  gas  exit  angle,  we  can  use  the  following 
empirical  equation  from  the  designed  black  profile  to  verify: 

Pecos'S-  (38) 

where  a  is  the  width  at  the  throat  of  the  outlet  of  the  duct 
(see  Figure  2).  The  obtained  result  is  P,  =  60.2*  .  It  is 
basically  consistent  with  the  design  requirement. 


Firure  ]  -  ■  *.  ♦  y  \ '  rt  ribution  on  the  blade  profile 

obtain-  :  .  lo  calculation. 

Key:  1.  ..  .u-.:,  r.i.-ar  fringe,  3.  calculated  point,  i( . 

real1  frli.g--,  •_>.  inner  arc,  6.  stagnation  point,  7.  back  arc, 

V.  The  Selection  of  Variation 
of  the  Design  Variables 

As  docori  :  -'U  above,  when  using  the  mean  stream  line 
method  to  design  planar  cascades,  the  design  variables  can 
be  chosen  :•  r ^  the  shape  of  the  mean  stream  line  and  the  velo¬ 
city  distribution  on  it  (which  includes  the  choices  of  the 
points  on  that  stream  line  that  intersect  with  the  front  and 
back  forehead  lines  and  by  how  much  the  velocities  at  those 
points  are  greater  than  the  given  velocity  values  in  the  front 
and  back  of  the  cascade).  Whether  a  good  cascade  can  be 
designed,  the  problem  lies  on  the  proper  selection  of  these 
design  variables. 

In  order  to  acquire  some  preliminary  experience  in 
designing,  we  have  carried  out  many  design  calculations  for 
various  types  of  blade  profiles  from  different  shapes  of  mean 
stream  lines  and  various  velocity  distributions  on  the  mean 
stream  line.  Because  the  solutions  to  the  reverse  problem  of 
the  subsonic  compressible  flow  is  qualitatively  very  close  to 
that  of  the  incompressible  flow,  therefore  during  the  calcula¬ 
tions  most  of  the  time  we  use  the  incompressible  flow  in  order 


to  obtain  more  information  for  comparison  more  conveniently. 

The  following  is  a  presentation  of  preliminary  recommendations 
on  the  selection  of  design  parameters  based  on  the  meaningful 
results  obtained  from  the  information  already  in  hand. 

1.  Pulsed  Cascade 

In  the  calculation,  a  comparison  is  carried  out  on  the 
different  velocity  distributions  of  pulsed  blade  family  with  a 
circular  mean  stream  line.  The  result  obtained  is:  Regard¬ 
less  of  whether  the  flow  is  compressible  on  incompressible  and 
the  distribution  of  velocity  on  the  mean  stream  line  (as  long 
as  the  variation  is  not  too  large),  the  trend  of  variation  of 
the  velocity  distribution  on  both  sides  of  the  duct  (that  is 
on  both  sides  of  the  blade)  is  similar  to  that  of  the  velocity 
distribution  on  the  mean  stream  line.  This  results  can  be 
qualitatively  explained  as:  because  the  curvature  of  a  circle 
is  uniform  everywhere,  therefore  the  difference  of  the  velo¬ 
city  at  the  boundaries  of  the  duct  is  not  too  different  from 
that  on  the  mean  stream  line.  The  typical  conditions  can  be 
seen  in  Figures  4  and  5*. 

From  these  two  figures,  it  is  easy  to  find  out  that  the 
aerodynamic  property  is  not  good  for  the  blade  profile  using 
a  circular  arc  as  the  entire  mean  stream  line.  Because  that 
the  felocity  distribution  of  the  back  arc  is  basically  parallel 
to  that  on  the  mean  stream  line,  in  the  segment  near  the  outlet 

*For  the  convenience  of  referring  to  the  blade,  the  velocity 
distribution  is  plotted  along  the  z  direction  rather  than 
along  the  surface  of  the  blade  profile  starting  from  Figure  4. 
Thus  the  absolute  value  of  the  gradient  of  velocity  appears 
larger  in  the  front  and  back  fringes  -  especially  at  the  rear 
fringe  of  the  back  arc.  The  actual  velocity  gradient  J  J 

Is  actually  much  small.  Please  notice  this  point  when 
referring  to  all  figures  afterwards. 


on  the  back  arc  of  the  blade  there  must  be  violent  compression 
and  gas  separation  caused  by  it.  In  the  internal  combustion 
turbine  laboratory  of  Tsung  hua  University,  we  have  carried 
out  a  lev;  velocity  winch  tunnel  experiment  on  the  planar  cas- 
case  basically  similar  to  the  design  shown  in  Figure  4.  The 
result  showed  that  there  is  a  separation  region  at  the  exit 
of  the  back  arc  of  the  blade.  The  flow  field  behind  the  cas¬ 
cade  is  extremely  inhomogeneous. 

When  placing  colored  liquid  in  the  separation  region  on 
the  back  arc,  we  can  use  our  naked  eyes  to  observe  the  appear¬ 
ance  of  whirlpools  in  this  region.  The  liquid  droplets  adhere 
to  the  blade  in  rotation  and  do  not  get  blown  away.  These 
indicate  that  the  design  using  circular  arc  as  the  entire 
mean  stream  line  is  not  good. 

But  for  ordinary  pulsed  cascades,  the  velocity  variation 
in  front  and  behind  the  cascade  is  not  large.  It  is  naturally 
proper  to  have  less  variation  in  the  velocity  distribution 
on  the  mean  stream  line.  There  is  not  too  much  choice.  The 
velocity  distribution  on  the  back  arc  at  this  time  is  usually 
desired  to  be  flat  and  invariant  at  the  beginning  and  then  it 
is  followed  by  a  segment  of  compression  region  with  a  velocity 
gradient  of  a  negative  value  which  is  not  too  large  and 
finally  connects  to  the  value  of  the  exit  velocity.  There¬ 
fore,  for  the  beginning  section  of  the  duct  in  a  pulsed  cas¬ 
cade,  the  use  of  a  circular  arc  as  the  mean  stream  line  and 
taking  the  value  of  the  velocity  as  a  constant  are  recommended. 
At  this  time  the  rising  section  of  the  back  arc  has  a  basically 
flat  distribution  of  velocity.  After  that  we  can  find  some 
curves  as  the  mean  stream  line  to  satisfy  the  latter  require¬ 
ment  discussed  above. 

In  order  to  compare  the  effect  of  the  shape  of  the  mean 


stream  line  more  conveniently,  the  blade  profiles  of  three 
types  of  mean  stream  lines  with  constant  velocity  distribution 
are  shown  in  Figures  C,  7,  ani  8.  The  first  halves  (fl/n>0’> 
of  the  mean  stream  lines  are  all  circular  arcs  and  the  rear 
halves  are  equi-angular  variation  curve  in  the  z 

direction  (its  analytical  equation  can  be  found  in  Table  2  of 
the  Appendix),  suspended  chain  curve,  and  parabolic  curve, 
respectively.  The  junction  point  of  various  mean  stream  lines 
is  fixed  at  o  because  only  here  it  is  convenient  to  make 
the  variations  of  (t,,  and  ft.  continuous. 

From  these  three  figures,  we  can  see  that  in  the  latter 
half  of  the  duct  where  ft,<  0  ,  the  effect  of  velocity  distri¬ 

bution  on  the  back  arc  for  the  three  mean  stream  lines  is 
different.  When  the  mean  stream  line  is  a  parabolic  curve, 
the  velocity  on  the  back  arc  initially  drops  rapidly  and  then 
gradually  flatens  out.  In  the  case  of  the  suspended  chain 
curve,  the  situation  is  the  uniformly  compression  along  the 
z  direction.  In  the  case  of  the  equi-angular  variation  curve 
In  the  z  direction,  the  compression  region  is  more  concentrated 
towards  the  rear  (but  as  discussed  above,  the  velocity  distri¬ 
bution  Is  plotted  in  the  z  direction.  When  transformed  into  a 
plot  along  the  surface  of  the  blade,  due  to  the  higher  devia¬ 
tion  away  from  the  axial  direction  in  the  rear  of  the  blade, 
the  negative  value  of  the  velocity  gradient  is  not  too  large). 

If  the  circular  arc  is  also  used  as  the  mean  stream  line  In 
the  ft„< o  region,  then  we  will  obtain  situations  in  which  the 
compression  region  is  concentrated  even  more  towards  the  rear. 
The  velocity  value  on  the  back  arc  is  basically  the  constant. 
Only  in  a  very  short  section  near  the  outlet  violent  compression 
occurs.  From  these  four  conditions,  it  seems  that  it  is  related 
to  the  variation  of  the  ft*  value  (the  corresponding  curvature) 
on  their  mean  stream  lines  after  ft.<  0  (refer  to  Table  2  in  the 
Appendix).  The  ft*  of  the  equi-z-directional  angular  variation 


curve  in  the  z  direction  is  a  constant.  The  absolute  value  of 
{£-- 0*  (v„-X)ot\v p„  for  a  suspended  chain  curve  is  becoming 
smaller  in  the  z-direction  (the  absolute  value  of  the  term 
)  is  getting  larger,  the  term  is  getting 

smaller,  but  the  latter  is  the  dominant  factor),  its  initial 
compression  this  arrives  faster.  The  absolute  value  of 
-  2  'I  cos:  p„,  for  the  parabolic  curve  becomes  smaller  even 

faster,  therefore,  the  compression  begins  even  more  rapidly. 

On  the  contrary,  the  absolute  value  of  — ~  •  for  a 

circular  arc  in  the  z  direction  is  getting  larger  when  pm< o 
The  onset  of  its  compression  then  is  slower  than  that  of  the 
equi-z-directional  angular  variation  curve.  -The  effect  of 
using  other  ordinary  curve  as  the  mean  stream  line  can  also 
be  evaluated  approximately  based  on  this  technique. 

From  the  above  figures,  it  can  be  found  that  for  pulsed 
cascades  we  can  recommend  a  circle  when  Pm<  0  and  an  equi-z- 
directional  angular  variation  curve  or  a  suspended  chain  curve 
when  o  as  the  mean  stream  line.  The  choice  between  the  two 

latter  ones  should  take  the  application  condition  of  the  cascade 
into  consideration.  If  the  estimation  shows  that  the  velocity 
distribution  of  the  former  will  not  cause  separation,  then  it 
is  better  to  use  the  equi-z-directional  angular  variation 
curve  because  the  density  can  be  reduced  for  the  same  cascade 
distance,  the  width  of  the  blade  can  also  be  decreased,  and  the 
frictional  loss  of  flow  can  be  minimized. 

As  discussed  above,  for  pulsed  cascades,  there  is  not  too 
much  choice  in  the  velocity  distribution  or  the  mean  stream 
line.  If  the  M  number  is  not  too  large,  we  can  directly  choose 
the  velocity  distribution  on  the  mean  stream  line  as  a  straight 
line  or  as  a  constant.  If  the  maximum  velocity  is  to  be 
reduced  and  the  compression  gradient  is  to  be  minimized  to  the 
extent  possible,  then  we  can  first  choose  the  lowest  possible 


Fig.  4.  The  blade  profile  and  velocity  distribution 
when  the  mean  stream  line  is  a  circle. 

(M.=  0,  0,=  SO*.  0,=  -6O\  — 

b  Oofti  * 

«*i-45,  —  -1.54,  — 3—  =  0.32,  -j- 
••0.47) 

Key:  1.  W  m/sec,  2.  center  line  of  the  inner  contact  circle 

of  the  duct,  3.  mean  stream  line. 

constant  as  the  velocity  on  the  mean  stream  line  and  then 
use  a  certain  curve  to  connect  to  the  value  at  the  rear 
fringe  (its  selection  will  be  shown  later).  For  example,  it 
has  been  done  this  way  in  the  above  example.  Another  example 
is  shown  in  Figure  9.  The  given  parameters,  the  shape  of  the 
mean  stream  line,  and  the  selection  of  the  cascade  distance 
of  this  blade  profile  are  the  same  as  those  in  Figure  7.  Only 
the  velocity  distribution  on  the  mean  stream  line  is  different. 
Comparing  these  two  examples,  we  can  see  the  effect  of  this 
change . 
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Fig.  5.  The  blade  profile  and  velocity  distribution 
when  the  mean  stream  line  is  a  circle. 

2.  P,  =  50\  0,=  -S5.2\  — “1.59,  -y-=1.5S,  — -O.J55,  “ 

“0.50) 

Key:  5.  W  m/sec,  6.  center  line  of  the  inner  contact 
circle  of  the  duct,  7.  mean  stream  line. 

2.  Reversed  Cascade 

There  is  little  experience  in  using  the  mean  stream  line 
method  to  design  the  reversed  blade  profiles^^.  Since  the 
mean  stream  line  and  the  velocity  distribution  on  it  are 
simultaneously  affecting  it  more  significant ly ,  it  is  then 
more  difficult  to  analyze  the  effects.  But  we  can  be  sure  that 
the  qualitative  effect  of  the  mean  stream  line  as  discussed 
above  is  still  approximately  cowrect.  For  the  comparison  of 
the  effect  of  velocity  distribution  on  the  mean  stream  line 
at  this  time.  Figures  10,  11,  and  12  show  the  results  of  three 
blade  profile  designs  with  different  velocity  distributions 
on  the  mean  stream  line.  (They  are  a  straight  line  H'm  =  tr0  +  , 

a  second  order  parabolic  curve  N',,  -  n.*  +  ij,r:  ,  and  an 

eights  order  polynomial  ir„  =if9-r  ’V-’i,**),  respectively),  but 
the  other  parameters  remained  the  same.  It  is  apparant 
that,  for  the  reversed  cascade,  the  velocity  gradient  on  the 
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Figure  6.  The  blade  pro-  Figure  7.  The  blade  profile  and 
file  and  velocity  dlstri-  velocity  distribution  when  the 
bution  when  the  mean  stream  mean  stream  line  is  a  circle  KP„>o) 
line  is  a  circle  and  a  suspended  chain  curve  (0„<o)  . 

and  an  equi^Zr-directional  Key:  5.  W  m/sec,  6.  center  line 

angular  variation  curve  of  the  Inner  contact  circle  of  the 

i(P»< °>  .  duct,  7  mean  stream  line. 

Key:  1.  W  m/sec, 

2.  center  line  of  the 
inner  contact  circle  of 
the  duct,  3.  mean  stream 
line . 


mean  stream  line  is  generally  positive.  From  these  three 
figures  and  some  computational  results,  we  found  that:  in 
order  to  make  the  compression  small  in  the  most  troublesome 


area  near  the  rear  fringe  of  the  back  arc,  we  should  choose 

the  smallest  W  '  value  at  that  place  possible.  W  "  is  a 
m  r  e  m 

negative  number  (Wm  may  be  relatively  larger  at  this  time) 


3/ 


as  shown  in  Pi  euros  3  and  12.  We  should  not  uc<-  the  smallest 

W  possible  in  front  of  the  rear  fringe  as  we  might  have 

imagined.  It  would  rapidly  become  larger  as  shown  in  Figure 

11.  Besides,  we  must  also  realise  that,  at  the  front  and  rear 

fringes,  in  addition  to  W  the  values  of  W  '  are  also  influ- 

m  m 

encing  the  thickness  of  the  front  and  rear  fringes  signifi¬ 
cantly.  The  front  (rear)  fringe  has  a  larger  W  '  than  it 
makes  the  from  (rear)  fringe  thicker.  This  can  be  found  by 
comparing  the  blade  profiles  in  Figures  10^12. 


Figure  8.  The  blade  profile  and  the  velooit-v  distribution  of  a 
mean  stream  line  which  is  a  circle  when  (3„>o)  and  a  parabolic 
curve  when  (B„<o)  •  Key:  1.  W  m/sec,  2.  center  line  of  the 
inner  contact  circle  of  the  duct,  3. mean  stream  line. 


Just  as  pointed  out  in  Reference  [33,  in  the  design  of 
reversed  blade  profile,  it  is  more  appropriate  to  use  the 
stream  line  which  has  60%  of  the  flow  away  from  the  back  arc 
as  the  mean  stream,  line.  When  the  analytical  method  is  used 
to  carry  out  the  design  work,  the  same  experience  has  been 
verified.  In  addition  to  the  fact  that  the  mean  stream  line 
is  farther  away  from  the  inner  arc  as  pointed  out  in  Reference 
[33,  this  is  also  due  to  the  better  convergence  towards  the 
back  arc  direction  when  using  Taylor  series  to  expand  Pir, 

Therefore,  in  the  design  of  blade  profiles  of  the  reversed 
type,  it  is  recommended  to  use  the  60£  stream  line  of  an 
equi-z-directional  angular  variation  curve  or  a  suspended 
chain  curve.  The  velocity  distribution  on  the  stream  line 
should  be  the  type  which  increases  uniformly  and  becomes  flat 
at  the  rear  fringe  for  better  results.  Another  example  of 
blade  profile  with  this  type  of  design  can  be  seen  in  Figure 
13. 


After  the  mean  stream  line  is  selected,  as  described 
above,  there  is  still  this  problem  left  to  determine  the  posi¬ 
tions  of  the  front  and  rear  fringes.  The  difficulty  to 
precisely  obtain  the  shapes  of  the  front  and  rear  fringes  to 
satisfy  the  requirements  of  entrance  and  exit  angles  of  the 
gas  is  one  disadvantage  of  the  mean  stream  line  method.  But 
based  on  the  experience  acquired  from  calculations,  it  should 
be  recommended  to  take  the  position  of  the  rear  fringe  at 

.  The  value  of  cos-*—  at  this  time  is  very  close  to 
that  of  6e-  For  all  the  calculations,  it  is  within  +  2°. 

If  the  gas  exit  angle  of  the  black  profile  must  be  changed 
by  In  order  to  satisfy  a  requirement  after  the  design 

work  is  completed,  then  it  is  recommended  to  preliminarily 
change  the  Bm  by  at  the  rear  fringe  as  an  approximation. 

Of  course,  we  can  also  find  the  exit  angle  of  the  gas  from  the 
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momentum  variation  in  the  duct  based  on  the  velocity  distri¬ 
bution.  But  because  the  velocity  distributions  at.  the  front 
and  the  rear  fringes  can  not  be  obtained  accurately  at  the 
present  time,  it  is  still  very  difficult  to  accurately  deter¬ 
mine  the  exit  angle  of  the  gas  using  this  method.  The  posi¬ 
tion  of  the  front  fringe  can  be  generally  taken  at 
This  is  more  appropriate  for  the  pulsed  cascades.  Because  of 

the  larger  gas  entrance  angle  8.,  at  that  time,  the  deviation 

1  T  4  1 

of  the  incoming  flow  by  the  cascade  is  smaller  .  In  addi¬ 
tion,  since  the  mean  stream  line  is  located  at  the  center  of 
the  duct,  the  circular  components  of  the  deviation  due  to 
blades  on  both  sides  can  easily  cancel  each  other.  Besides,  for 
pulsed  cascades  which  are  denser,  the  slight  change  in  the 
gas  entrance  angle  only  affects  the  area  near  the  front  fringe. 

The  effect  on  the  velocity  distribution  of  the  entire  blade 

r  c -| 

profile  is  very  small  .  Even  if  there  is  some  deviation,  it 
does  not  matter  that  much.  For  reversed  cascades,  the  position 
of  the  front  fringe  can  be  taken  as  t  ab  where  dp 

can  be  0°  ~  10°  preliminarily.  When  the  gas  entrance  angle  B. 

rm  1 

is  smaller,  then  dp  must  be  larger  . 

In  the  selection  of  the  velocity  distribution  on 

the  mean  stream  line,  the  W  values  at  the  front  and  back 

m 

fringes  are  major  factors  which  determine  the  thicknesses  at 
the  front  and  back  fringes.  For  conventional  blade  profiles, 
it  is  recommended  to  take  •  io~ i .  s.V  ii  ,  at  the  front  fringe 

and  ir„  =  (l.02~i.08)K'„  at  the  rear  fringe.  Naturally,  a 
large  value  of  Wm  should  be  used  at  the  front  or  rear  fringe 
when  thicker  fringes  are  desired.  In  addition,  the  effect  of 
Wmf  as  discussed  above  should  also  be  considered. 

In  the  design  of  the  cascades,  the  selection  of  cascade 
distance  can  be  made  by  referring  to  the  existing  cascade  data. 
But  if  the  distance  used  is  too  large,  it  may  not  be  possible 


Figure  9.  The  blade  pro-  Figure  10.  The  reversed  type 
file  and  velocity  distri-  blade  profile  and  velocity  distri¬ 
bution  of  the  mean  stream  bution  of  the  mean  stream  line 

line  which  is  a  circle  (0„>o;  which  is  an  equi-z-directional 
and  a  suspended  chain  curve  angular  variation  curve. 

Key:  1.  W  m/sec,  2.  Key :  5.  center  line  of  the  Inner 

center  line  of  the  inner  contact  circle  in  the  duct,  6. 

contact  circle  in  the  mean  stream  line,  7.  M  m/sec. 

duct,  3.  mean  stream  line. 

to  use  only  the  first  three  terms  of  the  Taylor  series  in  the 
mean  stream  line  method.  This  can  be  realized  by  noticing  a 
rapid  change  of  Pit'.,  in  the  y  direction  (usually  the  fact  that 
it  becomes  smaller  near  the  Inner  arc  shows  up  first)  when 
calculating  equations  (33)  and  (34).  Therefore,  in  the  present 
use  of  the  mean  stream  line  method  to  design  planar  turbine 
cascades,  it  is  more  appropriate  to  choose  a  cascade  density 
greater  than  1.4 

As  discussed  above,  the  solutions  of  the  blade  profiles 
of  the  subsonic  compressible  flow  and  incompressible  flow  are 


very  close.  With  the  sane  mean  stream  line,  basically  it  only 
shifts  the  blade  profile  towards  the  direction  of  the  inner 
arc.  (Strictly  sneaking  it  is  impossible  to  compare  between 
the  two  cases.  This  is  because  that  the  velocity  distribution 
on  the  mean  stream  line  cannot  be  the  same  for  the  compressible 
and  incompressible  cases  when  the  mean  stream,  lines  are  the 
same.  Similarly,  at  W^  the  compressible  W  is  larger,  there¬ 
fore,  the  velocity  distribution  values  of  the  blade  profile 
obtained  are  also  larger).  In  addition,  because  of  the  varia¬ 
tion  in  density,  when  the  W  of  the  front  (or  rear)  fringe  is 
greater  or  equal  to  that  of  W,  or  W  ,  the  increase  of  P«,  for 
the  compressible  case  is  smaller  which  means  that  the  front  (or 
rear)  fringe  can  be  made  thinner.  However,  this  effect  is  not 
significant.  Besides  these,  we  can  imagine  that:  in  the  sub¬ 
sonic  region,  it  may  be  more  accurate  to  use  the  mean  stream 
line  method  when  the  M  number  is  larger.  This  is  because  that 

the  trends  of  variation  of  W  and  p  are  opposite  where  using 

z 

Taylor  series  expansion  on  both  sides.  Thus  it  makes  it  poss¬ 
ible  that  the  series  of  PJf,  may  converge  more  rapidly. 

IV.  The  Design  Method  Which  Satisfies  a 
Given  Blade  Thickness  Distribution 

The  next  advantage  of  the  mean  stream  line  method  is  that 
it  allows  the  design  of  a  blade  profile  which  approximately 
satisfies  a  preset  thickness  distribution  in  order  to  ensure 

the  strength,  rigidity  and  cooling  requirements  of  the  blades 

rn 

.  When  using  tne  analytical  method  to  carry  out  the  cal¬ 
culation,  there  are  two  design  methods  to  assure  the  thick¬ 
ness  distribution  as  discussed  in  the  following. 

In  one  of  the  two  methods,  the  design  variables  are 
still  *  )  and  as  discussed  before.  In  order 

to  assure  the  distribution  of  thickness  at  this  time,  we 
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Figure  11.  The  reversed 
type  blade  profile  and 
velocity  distribution  when 
the  mean  stream  line  is  an 
equi-z-directional  angular 
variation  curve. 

Key:  1.  W  m/sec,  2.  center 

line  of  the  inner  contact 
circle  in  the  duct,  3.  mean 
stream  line  . 


Figure  12.  The  reversed  type 
blade  profile  and  velocity 
distribution  when  the  mean  stream 
line  is  an  equi-z-directional 
angular  variation  curve. 

Key:  5.  W  m/sec,  6.  60°/stream 

line . 


should  follow  the  instructions  in  Reference  [1]  to 
initially  assume  that  to  determine  the  approximate 

variational  curve  (PikJ  m~ z  .  After  the  shape  of  the  mean 
stream  line  '*«=*(*)  is  selected,  the  variation  curve  * 

can  be  obtained  from  the  first  order  isoentropic  aerodynamic 
equation.  That  means  we  can  use  an  analytical  expression 

*  )  which  is  very  similar  to  this  curve  to  replace  it 
to  solve  the  problem.  Since  the  relation  between  (Ptt'Ji  and 
•b-b  is  already  very  similar  and  the  thickness  distribution 
of  the  blade  in  actual  design  work  only  has  to  satisfy  the 
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=  1 . 79,  =  0.364*  O’ 

Figure  13.  The  reversed  type  blade  profile  and  velocity 
distribution  of  the  mean  stream  line  which  is  a  suspended 
chain  curve. 

Key:  8.  W  m/sec,  9.  60%  stream  line. 

requirements  approximately,  therefore  it  does  not  have  to  be 
too  rigorous  to  transform  (Pir..’‘.~  z  into  — ■  -  and  finally  onto 

w,im=ik(  3 )  in  order  to  avoid  mathematical  difficulties  and 
increase  the  computational  speed. 

The  other  method  also  used  the  selection  of  CP!Pj„=“(0 
and  y„=y(  z )  as  in  Reference  [1]  to  be  the  variables  in  order 
to  obtain  a  series  of  equations  from  which  the  analytical  solu¬ 
tion  can  be  found.  Their  forms  are  very  similar  to  the  ones 
given  in  Section  III  of  this  paper.  Since  in  the  actual  design 
it  does  not  require  the  strict  assurance  of  the  thickness 
distribution,  we  can  use  the  approximately  corresponding 

*  )  to  replace  (pk'«)«»=w(  z  )  as  the  variable.  The  actual 
equations  are  omitted  here  to  save  space.  It  has  been  docu¬ 
mented  in  detail  in  Reference  [8]. 
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VII.  Solving  the  Forward  Problem 
Using  the  Analytical  Method 

It  is  mere  complicated  to  solve  the  forward  problem  using 
the  analytical  method  because  it  is  necessary  to  find  the 
analytical  expressions  »',  =  »'(  z  )  and  y„,  =  y(  O  which  satisfy 

the  boundary  of  the  blade  profile  at  this  time.  However,  it 
still  has  some  advantages  over  the  numerical  differentiation 
method  used  in  the  past.  Since  it  is  unavoidable  to  carry  out 
a  few  trial  runs  to  approach  the  boundary  of  the  blade  profile 
in  solving  the  forward  problem^-  ^  ^  ^  ^ ,  it  is  of  course  more 

convenient  and  fast  to  perform  these  trial  runs  using  the 
analytical  method.  From  this  we  can  first  obtain  the  mean 
stream  line  and  the  trend  of  variation  of  parameters  on  it. 

We  then  can  use  the  numerical  solution  to  make  the  final 
correction.  Thus  the  entire  approaching  processes  can  be  com¬ 
pleted  faster  and  more  accurately.  Actually  for  those  who  are 
already  familiar  with  the  analytical  method,  it  is  more 
desirable  to  use  the  analytical  method  to  carry  out  the  entire 
approaching  process.  Based  on  the  experience  acquired  to  date. 
It  is  not  going  to  be  slower  than  using  the  numerical  method. 
Since  it  avoids  numerical  differentiations,  the  accuracy  is 
even  higher. 

Because  that  the  accuracy  in  using  the  analytical  method 
to  solve  the  reverse  problem  is  higher  and  there  are  lots  of 
calculated  blade  profiles  available  to  be  referred  to, 
therefore  we  have  verified  these  blade  profiles  in  order  to 
examine  the  validity  of  the  empirical  variational  rules  of 
*•  and  c  as  presented  in  References  [1]  and  [3]. 

The  result  of  verification  is  satisfactory.  We  can  only 
supplement  as  follows:  the  position  of  the  mean  stream  line 
for  an  incompressible  flow  Is  even  more  toward  the  inner  arc 
than  the  one  obtained  using  the  plotting,  method  in  Reference 
[3]  (the  physical  meaning  of  this  point  is  obvious).  We 
suggest  that  the  envelop  curve  **'b  in  Figure  5  of  Reference 
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[3]  bo  used  directly  as  the  m^an  strf'nn  lino*.  Tn  addition, 
under  many  conditions  there  are  similarities  between  the 

-  curve  and  the  icm$/p~~z  curve  (specific  situations 
are  shown  in  Reference  [8]).  Therefore,  in  solving  the  forward 
problem  using  the  analytical  method,  it  is  possible  to  use 
these  recommended  analytical  expressions  which  are  determined 
to  have  approximately  the  same  variations  as  y^  and  (pu',)Z 
or  (PilO-  to  obtain  the  preliminary  solution.  When  the  obtained 
duct  boundary  is  not  consistent  with  the  given  condition,  we 
can  carry  out  a  correction  as  presented  in  References  [1]  and 
[33-  and  then  list  the  analytical  expressions  after  the  correc¬ 
tion.  If  one  is  familiar  with  the  characteristics  of  various 
curves,  it  is  not  a  difficult  step  either.  As  discussed  above, 
in  the  selection  of  the  analytical  expressions  of  and 

or  it  is  possible  to  connect  them  with  several 
curves  which  is  more  convenient  in  finding  the  solution  to  the 
forward  problem.  If  it  is  difficult  to  make  H'*,,  ill,  Ill  and  s 
y  Pm.  Pm,  P™  continuous  at  the  junction,  it  is  also  possible 
to  separate  the  defined  region  by  assuming  that  another  segment 
of  a  smooth  curve  is  connecting  them. 

VIII.  Accuracy  of  the  Analytical  Method 

The  mean  stream  line  method  has  already  been  verified 

experimentally.  This  paper  merely  presents  an  improvement  in 

the  specific  calculation  procedure.  It  does  not  alter  its 

conditions  and  simplification  assumptions.  Therefore,  it  is 

not  necessary  to  conduct  another  experimental  verification. 

However,  for  the  further  verification  of  the  reliability,  the 

speed  of  solving  the  forward  problem,  and  the  applicability  to 

the  compressor  cascades,  we  have  carried  out  the  forward 

*  Comrade  We 1  Yu  Ping  who  is  a  graduate  student  of  Tsing  Hua 
University  Internal  Combustion  Machinery  Group  has  partici¬ 
pated  in  this  verification  work. 


problem  calculation  using  the  analytical  method  on  the  two 
cascades  with  experimental  data  shown  in  References  [4]  and 
[6],  The  comparison  of  calculated  and  experimental  results 
are  shown  in  Figures  14  and  15.  The  former  corresponds  to  a 
cascade  with  the  second  series  of  original  black  profile  in 
Reference  [A]  with  ;&  =  60*  t  attack  angle  (angle  between 
and  the  line  of  the  arc)  is  75°,  the  center  line  deflection 
angle  of  the  blade  profile  is  120°,  and  the  density  is  1.8. 

The  latter  is  the  Model  3  pulsed  compressor  planar  cascade 
in  Reference  [6].  The  purpose  of  this  calculation  is  to 
explore  the  feasibility  of  using  the  mean  stream  method  on 
an  axial  flow  compressor  planar  cascade.  References  [4]  and 
[6]  only  gave  the  experimental  data  of  the  above  two  cascades 
at  low  velocities.  Therefore,  in  our  calculation  we  assumed 
that  the  working  medium  is  an  incompressible  fluid. 

From  Figures  14  and  15,  it  can  be  found  that  the  calculated 
results  are  in  satisfactory  agreement  with  the  experimental 
data  (the  results  of  turbine  cascades  are  even  better.  For 
compressor  cascade,  although  the  absolute  values  of  the  velo¬ 
city  distribution  are  slightly  different,  yet  more  importantly 
the  trend  of  velocity  distribution  variation  is  completely 
consistent).  It  completely  satisfies  the  requirement  of 
engineering  calculation.  In  addition,  in  the  calculation  we 
have  verified  that  the  use  of  the  analytical  method  to  solve 
the  forward  problem  Is  effective  and  fast.  For  both  types  of 
cascades,  it  only  requires  about  five  trial  runs  to  fit  all 
the  coordinates  of  the  calculation  points  on  the  surface  of 
the  blade  profile.  Since  it  only  requires  very  small  amount  of 
time  to  complete  one  run,  therefore  the  total  required  time 
is  not  too  much. 


In  the  example  presented  in  Section  V  for  the  reverse 
problem,  when  the  mean  stream  line  is  a  circle  and  its  velocity 
distribution  is  a  constant,  the  flow  inside  the  duct  is  a  free 


Figure  1^.  Comparison  of  Figure  15.  Comparison  of  Computa- 
Computational  Example  of  tional  Example  of  the  Forward 
the  Forward  Problem  with  Problem  with  Experimental  Data. 
Experimental  Data.  Key:  4.  calculated  value. 

Key:  1.  calculated  5.  experimental  value. 

value,  2.  experimental 

value. 


whirlpool.  The  solution  at  this  time  can  be  obtained  rigor¬ 
ously:  The  boundary  of  the  duct  and  the  mean  stream  line  are 

concentric  circles  and  the  velocity  distributions  on  them  are 
also  constants.  Their  specific  numerical  values  can  also  be 
easily  obtained.  Comparing  the  results  of  the  rigorous 
solution  and  those  obtained  using  the  method  mentioned  in  this 
paper,  we  can  see  that  they  are  completely  consistent.  There 
is  only  very  small  difference  in  a  small  segment  of  the  front 
and  rear  fringes  and  those  areas  are  not  supposed  to  be 
accurately  determined  by  the  mean  stream  line  method.  This 
also  demonstrated  the  accuracy  of  this  method. 

Mi 


Finally,  it  should  bo  pointed  out  that:  the  above 
calculations  in  the  examples  although  were  carried  out  using 
ordinary  slide  rules  yet  the  coordinate  points  on  the  surface 
of  the  blade  profile  are  very  smooth  (when  slight  unsmoothness 
is  found  in  a  certain  location,  it  can  even  be  used  as  in 
indication  that  there  is  problem  in  calculation  at  this  point 
and  requires  an  examination).  From  this  the  advantage  of  the 
analytical  method  is  fully  demonstrated  (There  is  some  slight 
scatter  of  the  coordinates  obtained  using  the  numerical 
solution.  Refer  to  Reference  [3]). 

IX.  Conclusions  and  Projections 

This  paper  presented  an  improvement  procedure  for  the 
mean  stream  line  method  -  using  the  analytical  solution.  This 
improvement  can  significantly  simplify  and  speed  up  the  calcula 
tion  process.  It  does  not  require  the  multiple  usage  of  numeri 
cal  differentiations.  Slide  rules  can  be  used  in  the  calcula¬ 
tion.  The  work  load  required  to  correct  a  mistaken  term  is 
also  very  small.  It  also  simultaneously  raise  the  accuracy 
drastically.  In  addition,  we  can  control  the  design  variables 
to  calculate  good  blade  profiles.  Therefore,  this  method 
should  be  recommended  in  the  design  and  calculation  of  turbine 
cascades  of  usual  densities  and  pulsed  compressor  cascades  of 
high  densities. 

Other  than  that,  the  analytical  solution  itself  still  has 
a  lot  of  room  for  further  development  which  is  yet  to  be 
explored.  For  example,  using  its  high  calculation  speed 
characteristic,  we  can  obtain  a  series  of  cascades  for  compari¬ 
son  to  provide  detailed  information  on  the  selection  of  para¬ 
meters  in  the  design  of  cascades  or  to  even  directly  obtain  a 
family  of  cascades  to  be  selected  for  use.  Because  the  analy¬ 
tical  solution  is  used,  it  also  opened  the  possibility  of 
using  direct  analytical  method  (even  In  approximations)  to 


determine  tho  optimal  design  parameters.  When  it  is  necessary, 
the  analytical  solution  also  allows  the  conditions  that  higher 
terms  of  the  Taylor  series  expansion  with  respect  to  the  y- 
direction  must  be  used.  This  is  because  that  first  higher 
order  derivatives  can  be  more  precisely  determined  (it  is  not 
easy  to  obtain  accurate  values  using  numerical  differentia¬ 
tions)  and  secondly  the  use  of  more  terms  of  the  Taylor  series 
will  not  be  complicated  by  too  much  as  compared  to  the  use  of 
the  numerical  solution.  Thus  further  expansion  of  the  Taylor 
series  in  the  y-direction  can  be  used.  Then  it  is  possible  to 
design  general  compressor  cascades  and  to  select  the  inner 
(outer)  arc  of  the  turbine  cascades  to  calculate  the  outer 
(inner)  arc  using  the  mean  stream  line  method.  Especially,  it 
is  highly  probably  to  design  axial  compressor  cascades  with 
higher  density  near  the  inner  rdius  using  this  method.  Figure 
15  is  a  preliminary  example  of  such  a  case.  Although  the  den¬ 
sity  of  a  conventional  cascade  is  much  less  dense,  yet  its 
curvature  is  also  much  smaller.  This  corresponds  to  the  use 
of  a  small  segment  in  the  2  direction  of  the  blade  profile  as 
the  expansion  region  for  the  calculation.  Therefore,  it  seems 
that  it  is  even  possible  to  use  only  the  first  three  terms  of 
the  Taylor  series  to  calculate  compressor  cascades.  The 
capability  of  designing  one  blade  surface  from  the  surface  of 
one  blade  profile  of  the  turbine  cascade  also  provides  advan¬ 
tageous  conditions  to  assure  its  technological  and  aerodynamic 
properties.  In  addition,  it  is  also  worth  considering  as  to 
how  the  problems  on  the  arbitrary  revolving  surfaces  and  binary 
ducts  (such  as  the  axi-symmetrical  gas  inlet  tube  and  the  com¬ 
pressor  tube  of  turbine  cascades)  can  be  solved  using  the 
results  of  the  analytical  solution.  For  all  these,  we  wish 
that  development  will  be  gradually  make  in  the  future. 

Under  the  guidance  of  the  authors  and  Comrade  Chiang 
Shih  Yen,  students  of  the  1963  class  of  Tsing  Hua  University 
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gas  combustion  machinery  group,  Liu  Chin  Yu,  Shu  Weu  Tao, 
Show  Shau  Hau,  Wang  Chan,  Kuo  Kao  Tsian  and  Chen  Chun  Hua 
have  explored  part  of  the  problems  discussed  in  this  paper 
in  their  graduation  designs  to  various  extents.  Professor 
Wu  Chun  Hua  has  reviewed  this  paper  in  detail  and  often 
offered  valuable  suggestions  on  the  developmental  direction 
and  problems  in  this  work.  Comrade  Zhang  Chu  Chow  has 
plotted  all  the  figures  on  this  paper.  We  wish  to  express 
thanks  to  all  of  them. 
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APPENDIX 

The  following  two  tables  are  presented  in  the  appendix: 

1.  Table  2.  Table  of  Special  Equations  when  Frequently  Used 
Geometrical  Curves  are  Chosen  to  be  the  Mean  Stream  Line. 

2.  Table  3:  Table  of  Functions  and  Equations  of  Functions 
Used  in  the  Calculation. 


Table  2.  Table  of  Special  Equations  When  Frequently  Used 
Curves  are  Used  as  the  Mean  Stream  Line. 


Key:  2.  Colume  3-  Now;  4.  Geometric  Curve,  5.  Analyticaly 
relation,  6.  Circular  arc,  7.  Second  Order  Parabolic  Curve. 


Table  2  Cont’d. 


-  *  M  to  U  ? 


y  *=rp1  +  0i+  x 


'a/i-J(2n*  +  0) 


21)  -COS^fi 


—  8l)*cosa  (5  sin  fj 


Pl(2t\)U'  +  B1Hr> 


Pt (2  D  ) * •  W  +  /*i  ( J  tj ) \V>  +  B„ .  W' • 


P<-{2i\)t-(r  +  Ps(2r\)H''  +  Bll.IV+8u.'V'i 


w 


~trr 


-«jCB2.*(2H).IK  +  8s»^3 


(d^H'  \  r-  iy>2  T 

~dyi~)t  +  "sl/»<2  ’l  ):  W  +  P,(2  n  ) •*"  +  BltW  +  B^-^- J 
+  M,£(2  -ir)P»(2»))sH'  +  (7-3Y)-Pio(2’l)  W'-B-jv* 

W"a  n  wi’  t, 

+  (*3 1  -  Y  •*«)— Y  -*ao—  lK~ 


mm 


Table  2  Cent  rd 


Key:  1. 

Parabolic 
curve . 


x  / 
UA  ~ 

4 

! 

l 

~  dc  W  tt  «  sf 

2 

v  =v3+  X 

3 

4 

6n**cos:fJ 

1 

S  j 

t>»lco»s  f)  (  4  cos1  &  —  3  ) 

! 


7  r,(6D)»'  + 


8  i  *  •H’  +  W 


_ i _ w  _ 

10  1  (  *"-•  )  -.wc«so(6'i)*  •ic  +  fla«,'D 

|  \  oy  _ 

i  /  ^5|1  \  r  »"»  n 

I  V  ~a*»  +  +  »,  +  *l9--i7-J 

H  I  +  *{'  2  -  Y  >T«(1 2 1  > »'  »■  (7 -  3  Y *  >  «  •  *'  - 

I  1 

Wt  n  IV'2 

•  *(/<,,- Y»6)-~j-Mo.Y.BJ0 -y^,— 

.  a  >»•  i  r  /  ii',J  \a 

*  •  '  <#  -  ).  *MSL  ■•(«’•)*  >»"+«jj(»/''  +  ^-Y--77-)J 

Column,  2.  Now,  3-  (Continuation),  4  .  Third 
Curve,  5.  Equi  z  -  directional  angle  variati 


Table  2  Cont ' d . 


6 

1 

1  1 

&  a  « 

a  ! 

f 

y  ■=  T)ctw  h  (0^)+  X 

•  i 

tun_,C'lO»m  h  (0-:)] 

4  j 

e*(y  -  X)*->m3D 

5  j 

0-cos  0  sin  PC  1  ~  2  cos-  P  *0U(  y  —  X  )J1 

6 

HtQH>  -  X)l^ +  «;■•»" 

7 

i 

j  C0<(v  -  \)'-r?  +  QJ-nOn,  +  i‘i0i( y  -  + 

8  ' 

1  n.  <ijfv  _  x)if  +  »,<*" 

»'*> 


C<Wy  -  Us'\  +  a3-''d»''+/,i 0-(y  -  X)^' 

\v 


10 


(  d'Jy’  )  -X)w  +  W') 


(  ‘•j",-*-  )H  >  .  x  )j/>7_e«/j3iK+e»(y  -  \)r,w  +Biaw' 

+  1*19  ™u*-\+  w<[(  2  -  Y'On  »  -  x  )*./►„«'  + (7- JY)0»(V  -  X  y,o-K" 

n'»a  *1  i*"* 

*  (/  .I  -  V«b)  -K-J-A/6YWJ0-  — 

j j  !  ( od’r)fi + y  -  x >»" + «is( «' ♦«*•*  )] 


Key:  1.  Contlnuat ion ;  2.  Column,  3.  Row,  4.  Suspended  Chain 

Curve,  5>  Cosine  Curve. 


6° 


- - 


y  -  »lcos(8/)+  X 


— 

tnn^C  —  r)0sifi(O/)3 

-0J(  y  -  *)cos-0 

-0^x«psinpc  1  +  -  X)J3 

-PiO=(y  -  x)w  +  b.2w> 

C0‘(  y  -  A)2-P2-0-'/,]ty-Pj02(  y  -  A)H"  +  «„»'* 

-B13e=(y  -  A)if'  +  flai!" 

C04(v  -  A):P«-0-»:3ii'-/’i0-(y  -  x)W  +  Bi.,W'+Bn 

(-— -)B  +M=C«.<>e-(y  -  A)ir-B4iy"] 

)„  +  «i{c0*(  y  -  A)=P7  +  0V/j3iK-e3(  y  -  \)P,W  +  Bls»-» 

„  iytl  ^  r 

+  +  -X)l^-(7-SY)0 2(> 

ir»3  1  »"3 

-  /Wt"  +  (fl:  ( -  VH„)  — — J  -  M4Yft.0----  - 

fy 


Table  3.  Table 


)f  Functions  uO<_J  in  tin--  Calculation 

*3  rt  3 


B  ft  cosf>  *•"  0 


,  cosP(2cs»:p-  1)  2<x»:  pain  p 


sfttta  -i  I  * 


cos  P  (12cos«  p;  cos=  p  Sin  0  2  cos:  p  sin  p 

j— ncos>p  +  l)]*(4cos:p-i)  x  (geos5  P- J) 


1.CC00 

o.ocoo 

.0(300 

0.0000 

2.CC00 

0.9962 

0.08716 

0.9811 

0.1730 

1.8947 

0.9348 

0.1736 

0.9254 

\ 

0.3368 

1.5943 

0.9659 

0..dS8 

0.8365 

0.4830 

1 

1.1427 

0.9397 

0.3420 

• 

0.7198 

0.6040 

0.6047 

0.9063 

0.4226 

i 

0.5826 

0.6943 

0.05521 

0.6660 

1 

1  0.500C 

0.4330 

0.7500 

-0.4330 

0.8192 

0.5736 

| 

0.2802 

0.7698 

-0.8012 

0.7660 

0.6423 

0.1330 

0.7544 

-1.0133 

0.7071 

0.7071 

0.0000 

0.7071 

-1.0607 

0.6428 

0.7660 

-0.1116 

0.6330 

-0.9618 

0.5736 

0.8192 

-0.1962 

0.5390 

-0.7572 

0.500C 

0.8660 

-0.2500 

0.4330 

-0.5000 

0.4226 

0.9063 

-0.2717 

0.3237 

-0.2459 

0.3420 

0.9397 

-0.2620 

0.2198 

-0.04192 

0.2588 

0.9659 

-0.2241 

0.1294 

0.08204 

9.1736 

0.9818 

-0.1632 

0.05939 

0.1179 

[SBbi 


3?3  rt  n  m  a 


Efttt  5 


0.07548 

0.04001 

0.01734 

0.005235 


0.0000 

0.08584 

0.1533 

0.2253 

0.2667 

0.2851 


0.2046 

0.2214 

0.1768 

0.1768 

0.3364 

0.1308 

0.4493 

0.08866 

0.4063 

0.05413 

■ 

0.3046 

0.02891 

0.1686 

0.01286 

0.09125 

0.001334 

0.02368 

0.0008951 

0.0000 

-1.7128 

-3.2362 

-4.4091 

-5.1217 

-5.3307 

-5.0625 

-4.4054 

-3.4919 

-2.4749 

-1.5011 

-0.6882 

-0.1083 

0.2195 

0.3252 

0.2755 

0.1549 


coa>  0 


s' 

« 


1.0000 

0.9811 

0.9263 

0.8409 

0.7327 

0.6115 

0.4871 

0.3670 

0.263* 

0.1768 

0.1097 

0.06208 

0.03125 

0.01348 

0.004680 

0.001161 

0.000157# 


S5E 


is.  a  (-hi) 


*1! 

Cl 

c4  1 

C* 

e» 

cosJ0  x 

(«-7cos>0) 

1-2  cos  J  0 

i  ! 

i 

1  COS  psi  np  k 

I(4cos3  0  -3)i(3-  8COS-  0)! 

3  cm-  -  1 

1 

sin  0  x 

(1-6CO6*0) 

cos  0  x 

(3-acosJ  0) 

cos  0  sin  0  x 

(26cosi0-6) 

1  !  1 

n 

< 

a 

*  i 

a  | 

£  | 

*  i 

1 

-S' 

a 

£ 

0 

r 

a 

«f 

0.09590 

0.03011 


2.0000 

1.9319 

1.7321 

1 . 4142 

1.0000 

0.5176 

0.0000 

-0.517S 

-1.0000 

-1.4142 

-1.7321 

-1.9319 

-2.0000 

-1.0319 

-1.7321 

-1.4142 

-1.0000 


0.0000 

-0.8577 

-1.6278 

-2.2321 

-2.6124 

-2.7357 

-.25981 

-2.2233 

-1.6688 

-1.0000 

-0.3008 

0.3459 

0.8660 

1.2038 

1.3268 

1.2321 

0.9136 


0.7605 

-1.6204 

0.3003 

-1.4142 

0.2393 

-1.1330 

-0.01303 

-0.7978 

-  0.25’0 

1 

-  0.4330 

I 

-  0.1612 

i  -0.06492 

-0.6191 

) 

0.2802 

-0.7990 

|  0.5777 

1 

-  0.9095 

0  8066 

-5.0000 

-4.9204 

-4.6855 

-4.3120 

-3.8191 

-3.2366 

-2.5981 

-1.9398 

-1.2981 

-0.7071 

-0.1963 

0.2111 

0.5000 

0.6610 

0.7060 

0.6378 

0.4791 


0.0000 

1.7193 

3.2861 

4.5646 

5.4504 

5.8818 

5.8457 

5.3780 

4.5384 

3.50C0 

2.3353 

1.1999 

0.2165 

-  0.5195 
-0.9509 

-  1.0646 

-  0.3920 


5.  Even,  6.  Odd 


m  ft 

•  1 

ft 

ft 

ft 

ft 

f% 

ft 

&  * 

COS*  3  *  1 

cos*  f)  (20cos‘  (1 

2  cos*  0  sin  (}  x 

cm*  0  x 

cos*  9  sin  &  * 

cos*  (JsiajJ 

i 

(2cosip  -  J) 

-21cos!fl  (•  3) 

C  3  cos>  fj  -  3  ) 

(  9  cos*  0  -  7  ) 

(  8  cos*  0  -  l  ) 

; 

4+ 

/ 

f 

& 

m  i 

* 

" 

« 

fr-o*  1 

1.0000 

2.0000 

0.0000 

2.0000 

0.0000 

0.0000 

S* 

0.9736 

1.8218 

0.8479 

1.8880 

0.4269 

•  •  0.08551 

0.8975 

1.3388 

1.5S45 

1.5772 

0.7993 

0.1609 

15* 

0.7805 

0.6870 

2.0116 

1.1347 

1.0725 

0.2176 

20* 

0.63S6 

0.03745 

2.1677 

0.6522 

1.2198 

0.2506 

IS* 

0.4785 

-0.1620 

2.0365 

0.2400 

1.2359 

0.2584 

SO* 

0.32-18 

-0.73C7 

1.6875 

-0.1055 

1.1367 

0.2436 

3S* 

0.1871 

-0.7657 

1.2231 

• 

-0.2903 

0.9294 

0.2105 

40* 

0.07806 

-0.6426 

0.7302 

-0.3473 

0.7284 

0.1696 

45* 

0.0000 

-0.4119 

0.3536 

-0.3125 

0.5000 

0.1250 

50* 

-0.04612 

-0.2483 

0.07989 

-0.2315 

0.3009 

0.08106 

55* 

■^-0.06454 

-0.1083 

-0.06527 

-0.1438 

0.1505 

0.05085 

•0* 

-0.06250 

-0.03125 

-0.1083 

-0.07122 

0.05413 

0.02706 

•5* 

-0.04852 

-0.001520 

-0.09085 

-0.03072 

0.001901 

0.01222 

70* 

-0.03065 

0.003824 

-0.05309 

-  0.009520 

-0.01120 

0.004398 

75* 

-0.01502 

0.00(955 

-0.02138 

-0.001923 

-0.01002 

0.001218 

80* 

-0.001920 

0.0003768 

-0.001841 

-0.0001845 

-0.004224 

0.0001555 

Key: 


1.  Function,  2.  Equation,  3.  Functional  Characteristics , 
4.  Even,  5.  Odd. 
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-5.0000 

0.0000 

1 

1.0000 

0.0000 

0.0000 

; 

0.0000 

-4.7426 

-1.6938 

0.9663 

0.08519 

0.5737 

0.5014 

-4.0314 

-3.1386 

0.8713 

0.1584 

0.9424 

0.8778 

-3.0292 

-4.1137 

0.7320 

0.2102 

0.9837 

1.0434 

-1.9484 

-4.5226 

0.5713 

0.2355 

0.7025 

0.9759 

-0.9863 

-4.3786 

0.4126 

0.2342 

0.2209 

0.7197 

—  0.2740 

-3.7969 

0.2740 

0.2109 

-0.2813 

0.3654 

0.1460 

-2.9561 

0.1653 

0.1733 

-0.6426 

0.01646 

0.3133 

-2.0491 

0.09084 

0.1299 

-0.7802 

-0.2437 

0.3094 

-1.2374 

0.04419 

0.08839 

-0.7071 

-0.3750 

0.2256 

-0.6202 

0.01873 

0.05403 

-0.5064 

-0.3836 

0.1309 

-0.2264 

0.006719 

0.02917 

-0.2813 

-0.3091 

0.06055 

-0.02706 

0.001953 

0.01353 

-0.1083 

-0.2030 

0.02159 

0.03921 

0.0003301 

0.005164 

-0.01477 

—  0.1073 

0.005481 

0.03804 

0.00006404 

0.001504 

0.01417 

-0.04351 

0.0008450 

0.01846 

0.000005212 

0.0002904 

0.01142 

-0.01211 

0.00005469 

0.004671 

0.0000001436 

I 

0.00002700 

1 

0.003484 

-0.001782 
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Table  3«  Tabic  of  Functions  Used  in  the  Calculations  (concentration 
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function  r, 

i 

T , 

»» 

equation 

cos*  0  sio  0  x  ! 

(21  —  34COs.:  0 ) 

> 

cos*  0  sin  0 

>  cos>  0  sin2  0  x 

t 

|  (4  cos)  0-1) 

|  2  cos*  0  sin*  0 

3  cos*  0  sin>  0 

functional 

characteristics 

odd  i 

odd 

l 

even 

i 

even 

even 

=  0* 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

s* 

-1.0854 

0.08454 

0.02230 

0.01496 

0.02236 

10* 

—  1.8969 

0.1336 

0.08293 

0.05673 

0.08360 

15* 

-2.2340 

0.1961 

0.1649 

0.1166 

0.1690 

20* 

-2.1218 

0.2079 

0.2(38 

0.1824 

0.2571 

25*  J 

-1.6225 

0.1304 

0.3039 

0.2410 

0.3276 

30* 

-0.3492 

0.1582 

0.3248 

0.2813 

0.3654 

35* 

-0. 2237 

0.1163 

0.3045 

0.2963 

0.3623 

40* 

0.1361 

0.07623 

'  0.2502 

0.2846 

0.327Q 

45* 

0.3536 

0.04416 

0.1768 

0.2500 

0.2652 

50* 

0.3756 

0.02232 

0.1017 

0.2004 

0.1932 

55* 

0.2863 

0.009596 

0.04001 

0.1453 

0.1250 

60* 

0.1692 

0.003383 

0.0000 

0.09375 

0.07031 

65* 

0.07708 

0.0009223 

-0.01771 

0.05241 

0.03322 

70* 

0.02361 

0.0001760 

-0.01880 

0.02417 

0.01240 

75* 

0.005436 

0.00001945 

-0.01184 

O.OC8373 

0.003251 

60* 

0.0005333 

0.0000008142 

-0.004466 

0.301764 

0.0004S94 
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TIIE  ANALYTICAL  SOLUTION  OK  MEAN-STREAM  LINE 
METHOD  FOR  TWO-DIMENSIONAL  CASCADES 

- SOME  DEVELOPMENT  OF  MEAN- 

STREAM  LINE  METHOD C I ) - 

Tsay  Ruey-shen 
ABSTRACT 

A  modification  of  mo-.n-strcum  line  method  is  presented.  With  this  modification*  ait 
computation  of  two-dimensional  cascade  can  be  calculated  analytically  and  the  accuracy 
solution  is  higher  than  the  original  method  as  well  as  the  time  needed  for  one  complete 
calculation  is  shorter,  home  examples  of  solution  are  compared  and  discussed  for  finding 
the  better  form  of  the  mean-stream  line  and  the  bette-  variation  of  velocity  alone  the 
mean-stream  line.  With  this  discussion  and  the  equations  and  tables  presented  in  this 
paper*  only  abrut  S  nours  arc  needed  fir  projecting  a  good'  turbine  cascade  circulated 
by  compressible  flow,  and  the  Iteration  in  the  direct  problem  to  fit  the  Riven  blade 
shape  is  more  practical.  "1  he  results  of  calculation  of  this  method  compared  well  with 
experimental  data.  Some  probabilities  of  development  of  the  mean-stream  line  method 
are  discussed  in  the  end  of  this  paper. 
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CALCULATION  OF  TRANSONIC  FLOW  IN  PLANAR  CASCADES 
SOME  DEVELOPMENT  OF  MEAN  STREAM  LINE  METHOD  (IV) 


Tsay  Ruey-shen 


Abstract 


Because  the  equations  of  the  mean 
stream  line  method  are  consistent  at  sub¬ 
sonic  and  supersonic  speed  and  there  does 
not  exist  an  odd  point  at  M=l,  therefore 
they  can  be  used  to  calculate  transonic 
flow  in  the  cascades.  Comparing  the  cal¬ 
culated  results  with  the  experimental  data, 
we  found  that  it  is  feasible  to  use  this 
method.  The  approximate  range  of  condi¬ 
tions  under  which  this  method  is  applicable 
is  also  analysed  in  this  paper. 
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Symbols 


sonic  speed 
function  of  8l 
coordinate  in  the  direction  of  the  stream  line 
Mach  number  of  the  gas  flow 
direction  normal  to  that  of  the  stream  line 
distance  of  expansion  in  the  n  direction 

relative  flow  density  of  the  gas  flow 

<J*  i 

velocity  of  the  gas  flow 

coordinate  tangent  to  the  direction  of  the  cascade 
distance  of  expansion  in  the  y  direction 
axial  coordinate  of  the  turbine  machinery 
arc  tg  ,  gas  flow  angle 


ts<a£> 


ratio  of  specific  heats  of  gases 
W/a,  relative  velocity  of  the  gas  flow 
density  of  the  gas 
flow  function 

amount  of  flow  within  the  distance  of  expansion 


( 


fl  radius  of  curvature  of  the  stream  line 

Superscript 

d/dz,  total  differentiation  with  respect  to  the  z 

direction 

Subscript 

component  in  the  1  direction 
value  on  the  mean  stream  line  (selected) 
value  on  the  inner  arc  of  the  blade  profile 
value  on  the  back  arc  of  the  blade  profile 
component  in  the  z  direction 
value  under  one-dimensional  assumption 
n  different  functions 

critical  values  of  the  gas  flow 

I.  Introduction 

In  recent  years  transonic  cascades  have  been  widely  used 
in  turbine  machinery,  it  is  there  imperative  to  obtain  better 
understanding  about  the  state  of  flow  in  those  cascades.  It 
is  well  known  that  even  under  the  simplified  condition  of  an 
ideal  gas  flow  in  a  planar  cascade,  there  is  a  certain  degree 
of  mathematical  difficulty  in  solving  for  the  transonic  flow. 
This  is  because  that  the  basic  differential  equations  des¬ 
cribing  the  flow  in  the  subsonic  and  supersonic  regions  are 
different  types.  One  is  the  elliptical  type  while  the  other 
is  the  hyperbolic  type.  It  is  not  convenient  to  treat  them 
in  an  unified  manner. 

But  the  above  problem  only  appears  in  multi-dimensional 
equations.  In  one-dimensional  flow,  the  subsonic  and  super¬ 
sonic  flows  are  both  represented  by  an  integrated  algebraic 
equation  with  the  exception  that  we  must  notice  that  the  M 
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number  in  a  double  value  function  of  the  flow  density.  There¬ 
fore,  we  can  assume  the  following:  The  difficulty  in  solving 
the  transonic  flow  problem  is  because  the  multi-dirr.ensional 
flow  equation  is  used  directly.  If  we  can  initiate  the  general 
potential  flow  theory^^  of  turbine  machinery  to  continuously 
decrease  the  dimensions  all  the  way  to  one  then  we  can  use 
the  simple  method  to  directly  obtain  the  solution  of  transonic 

r  2 1 

flow.  The  mean  stream  line  method  J  is  sucb  a  method  in 
nature.  It  transforms  the  two-dimensional  flow  into  a  flow 
on  the  one  dimensional  flow  on  the  mean  stream  line  y  -y  (z) 
and  the  expansion  of  a  flow  parameter  in  another  direction 
(usually  in  the  tangential  direction  y  of  the  turbomachine. 

Thus  it  may  be  possible  to  solve  the  transonic  flow  problem 
using  an  unified  mathematical  treatment.  Especially  if  the 
analytical  solution  method  is  used  instead  of  the  numerical 
method,  it  is  much  more  easier  to  verify  this  possibility  and 
more  suitable  for  the  kind  of  calculation  accuracy  required 
in  the  transonic  region. 

Based  on  the  above  consideration,  this  paper  presents  a  dis¬ 
cussion  of  solving  the  transonic  flow  in  planar  cascades  using 
the  mean  stream  line  method.  Detailed  introduction  of  this 
method  has  already  been  documented  in  References  [2]  and  [3], 
so  it  is  not  repeated  here. 

II.  Equations  of  Mean  Stream  Line  Method 
Applicable  to  Transonic  Flow 

Reference  [3]  has  derived  the  equations  of  the  mean 
stream  line  method  which  are  suitable  for  the  analytical  solu¬ 
tion.  But  in  order  to  be  applicable  to  high  speed  flow,  it 
is  better  to  alter  the  forms  slightly.  In  addition,  in  order 
to  define  the  boundary  of  the  flow,  we  still  must  solve  a 
third  order  algebraic  equation.  Finally,  we  can  use  the  flow 
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function  concept  to  derive  a  form  which  is  even  easier  to 
solve  (see  Appendix  A).  It  is  summarized  in  the  following: 


The  partial  derivative  expanded  in  the  direction  tangent 
to  the  mean  stream  line: 
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The  equation  which  determines  the  flow  boundary  is 
(notice  that  Aip  has  plus  and  minus  signs): 
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The  equation  which  determines  the  velocity  distribution 
is 
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In  the  practical  calculation,  usually  the  remaining  small 
terms  in  equations  (5)  and  (6)  are  neglected.  The  B  '  s  in 
equations  (1)  -  (4)  are  simple  functions  of  6  and  their  defi¬ 
nitions  and  numerical  values  can  be  found  in  Reference  [3]  or 

in  Appendix  B.  Furthermore,  8*  and  6”  can  be  obtained  from 

’mm 

W  ( z )  as: 
m 


III.  The  Feasibility  of  Solving  Transonic  Flow 
Without  Intense  Shock  Wave  Using  the  Mean 
Stream  Line  Method 

As  discussed  above,  the  mean  stream  line  method  reduces 
a  two  dimensional  planar  problem  to  a  one  dimensional  flow  on 
the  mean  stream  line  and  the  expansion  of  a  flow  parameter  in 
another  direction.  There  is  not  difficulty  and  doubt  in  cal¬ 
culating  the  one  dimensional  transonic  flow.  Therefore,  the 
only  thing  which  we  must  check  in  solving  for  the  transonic 
flow  is  that  whether  it  is  feasible  to  expand  into  a  Taylor 
series  in  another  direction  which  intersects  the  stream  line. 

Apparently,  if  at  places  where  there  is  obvious  jumping 
variations  in  the  flow  parameters  it  is  impossible  to  obtain 
the  solution  using  the  Taylor  series  expansion.  Therefore, 
the  mean  stream  line  method  cannot  be  used  in  regions  with 
strong  shock  waves.  Besides  this  point,  because  the  same  type 
of  basic  equations  are  used  in  both  the  subsonic  and  supersonic 
regions  in  the  derivation  of  these  partial  differential  equa¬ 
tions  used  in  the  expansion,  the  equations  obtained  have  no 
odd  point  at  M  =  1,  as  long  as  the  flow  field  varies  relatively 
uniformly  and  smoothly,  it  is  as  effective  for  subsonic  flow 
as  for  supersonic  flow  when  using  the  series  expansion  method 
to  determine  the  corresponding  flow.  Since  it  is  possible  to 
use  the  mean  stream  line  method  for  the  form,  the  same  method 
should  also  be  applicable  to  the  latter.  For  good  cascades 
the  flow  inside  basically  varies  relatively  uniformly  and 
smoothly  which  coincides  with  the  above  requirement.  The 
exceptions  are  areas  near  the  front  and  rear  stationary  points 
where  the  mean  stream  line  method  cannot  be  used  to  obtain 
accurate  solutions. 


In  addition,  we  may  consider  the  following:  The  potential 
flow  is  a  reversible  flow  and  the  boundary  condition  of  planar 
cascades  is  usually  that  homogeneous  flow  exists  at  infinite 
distance  both  upstream  and  downstream.  Therefore,  when  it  is 
reduced  to  a  one-dimensional  problem,  the  relation  between 
each  flow  parameters  in  the  direction  of  the  expansion  is  also 
a  one-dimensional  flow  relation  (passing  through  one  stream 
line  to  infinity  and  then  flowing  back  on  another  stream  line). 
There  is  no  difficulty  involved  such  as  the  ones  encountered 
in  the  treatment  of  multi-dimensional  transonic  problems.  The 
only  thing  is  that  the  flow  variables  (there  is  only  one  var¬ 
iable  in  the  one  dimensional  isoentropic  flow)  are  not  defined 
through  the  use  of  variations  of  parameters  in  the  direction  of 
the  mean  stream  line. 

In  the  one-dimensional  flow,  the  only  difference  between 
the  supersonic  flow  and  the  subsonic  flow  is  that  the  flow 
density  of  the  former  is  a  decreasing  function  of  M  number  and 
latter  is  an  increasing  function.  When  M  =  1,  q  =  1  reaches 
the  maximum.  Therefore,  in  the  expansion  to  obtain  the  solu¬ 
tion,  It  is  only  natural  to  consider  that  these  relations 
should  be  checked  to  see  whether  thev  are  satisfied.  In  prac- 
tice,  when  M  =  1,  we  should  get  ii  *—  q  and  Jr'.-i  <n  .  When 

M  <  1,  *  should  have  the  same  sign  as  .  When  M  >  1, 

v'1  on, 

they  should  have  opposite  signs.  Using  equations  (A  -  10)  - 

(A  -  12)  in  Appendix,  we  can  see  that  these  relations  are  all 
satisfied . 

T  41 

There  w as  a  successful  example1-  J  using  ries  expansion 
in  a  direction  which  intersects  the  stream  line  to  solve  for 
the  transonic  straight  axial  nozzle  which  indicates  that  the 
above  idea  is  feasible. 


Of  course,  there  is  another  problem.  In  the  classical 


subsonic  mean  stream  line  method  it  is  only  necessary  to  use 
up  to  second  order  partial  derivative  of  velocity  in  the  Taylor 
series  expansion.  Is  it  enough  to  use  up  to  the  second  order 
term  in  the  transonic  region?  In  other  words,  when  only  up  to 
the  second  order  derivatives  are  used,  how  far  can  the  expan¬ 
sion  distance  be?  Can  this  be  applicable  to  general  transonic 
cascades?  This  problem  will  be  examined  through  the  actual 
calculation  and  qualitative  analysis  presented  in  the  next 
section. 


IV.  Expansion  Terms  and  Effective 
Expansion  Distance 

The  mean  stream  line  method  is  a  simple  method  which 
enables  us  to  obtain  results  in  short  period  of  time  using 
slide  rules.  Therefore,  although  in  principle  the  more  terms 
we  use  the  more  accurate  it  is,  yet  in  the  actual  application 
in  subsonic  cascades  we  only  used  up  to  the  second  order  par¬ 
tial  derivatives  of  velocity.  Of  course,  we  wish  that  we  can  do 
the  same  thing  for  the  transonic  flow.  The  calculated  results 
indicate  that  it  is  possible  for  ordinary  cascades.  The  examples 
of  the  calculation  process  are  shown  in  Figures  1  and  2. 

Figure  1  is  the  NACA  turbine  cascade  whose  experiment  results 
are  obtained  from  Figure  42(b)  in  Reference  [5].  In  Reference 
[6],  it  has  carried  out  a  calculation  for  this  cascade  using 
a  time  related  method.  The  results  are  also  plotted  here. 

From  this  figure,  it  can  be  found  that  the  hand  calculation 
obtained  using  the  mean  stream  line  method  is  satisfactory. 

It  seems  to  be  even  better  than  the  results  obtained  using  the 
time  related  method  by  an  electronic  computer  for  this  speci¬ 
fic  example.  Figure  2  shows  the  calculated  result  for  the 
TP-1A  turbine  cascade  of  which  is  also  satisfactory. 

But  trial  calculation  indicates  that,  when  the  Mm  number 
In  actual  working  condition  is  larger  and  A^  is  also  larger. 
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of  velocity.  Or,  in  other*  worix,  ;•  i :  1  •  ••s  that.  its  effec¬ 
tive  expansion  distance  is  ; — ;  ♦  *  •  r  »•••*.  This  pheno¬ 

menon  can  be  explained  by  t  !.„•  :.atu:v  tf  the  an  stream  line 
method  anJ  the  expansion  o.juut  i^:.s  ( I'.uuat  1  v.  :.x  { -12)  -  (A-15) 
in  the  Appendix)  in  the  normal  direction  of  it  s  velocity. 
Because  this  method  is  the  sim.r  lificaticn  of  a  two-dimensional 
flow  to  a  one-dimensional  flow  in  the  direction  of  the  mean 
stream  line  plus  the  correction  in  another  direction  (which  is 
reflected  by  the  partial  derivative  of  the  velocity  in  that 
direction).  Therefore,  the  more  the  flow  approaches  a  one- 


Figure  1.  The  cascade  pattern  and  its  experimental  and 
calculated  results.  Key:  1.  relative  width  of  the  cascade, 

2.  experimental  data,  3.  mean  stream  line  method,  4 .  time 
related  method. 

dimensional  flow  or  the  smaller  the  value  of  the  partial 
derivative  (the  less  correction  required),  the  more  accurate  it 
is;  which  also  means  the  farther  the  allowable  expansion  dis¬ 
tance  is.  Of  course,  when  the  flow  curvature  _1_  with 

hi 


respect  to  the  width  of  the  flow  boundary  (corresponding  to 
the  width  of  expansion)  and  the  expansion  and  contraction  of 
the  flow  boundary  ( correspond  ins  to  ('1'^  )  are  smaller,  the 


Figure  2.  Cascade  pattern  and  its  experimental  and 
calculated  results.  Key:  6.  relative  width  of  the  cascade, 
7.  experimental  data,  8.  mean  stream  line  method. 


more  close  to  one-dimensional  the  flow  becomes.  Therefore, 


the  application  of  the  mean  stream  line  method  should  be 

AV'  ,  A'X/ 

tor 


limited  to  the  situation  that 


t  S  \  • 

">r  i;  y\ 

A  v*  v  .... 


(or  )  cannot  be  too  large.  This  point 

is  also  shown  in  the  expansion  equations  (A-12)  -  (A-15) 
because  the  velocity  correction  due  to  the  presence  of  some 
partial  derivatives  consists  of  terms  with  and  '-f —*}  .in, 

(the  flow  parameter  variations  of  a  good  cascade  is  smooth 
so  that  only  lower  order  derivates  are  considered.  From 
the  expression  of  the  first  order  derivative  (A-12),  it  can 


)  and 


found  that  it  is  necessary  to  require  that  the  term  ,/,o 

which  is  affecting  the  curvature  cannot  be  too  large.  If  we 
use  this  method  to  obtain  the  solution  of  a  free  whirlpool 
with  variation  of  direction  but  no  variation  in  speed  and 
compare  the  result  with  the  rigorous  resolution,  we  realise 
that  in  order  to  obtain  the  necessary  engine^rl qcr  accuracy 
usually  it  is  desired  to  have  —  <  lit  or  --- —  $0  3z> 

(the  ordinary  good  cascades  are  aii  within  tms  range).  There 
is  another  point  worthwhile  our  attention  which  is  that  the 
limitation  in  the  magnitude  of  "f1  is  not  related  to  the 

Mm  number  of  the  flow  (see  Equations  (A-12)  and  (A-13)). 

Indeed,  the  velocity  distribution  of  a  free  whirlpool  is  not 

related  to  the  M  number.  Therefore,  when  we  use  the  mean 

stream  line  method  to  solve  the  free  whirlpool  or  a  flow  which 

principally  involves  a  change  of  direction  (corresponding 

impulsive  cascades).  Mm  has  no  significant  effect  on  the 

accuracy  of  the  calculation  which  means  that  there  is  no  added 

difficulty  in  solving  for  the  transonic  flow.  However,  on 

the  contrary,  from  equation  (A-13)  it  can  be  found  that  the 

effect  of  the  acceleration  ,  on  the  mean  stream  line  on  the 

«u  ». 

accuracy  of  the  calculation  is  closely  related  to  Mm  and  it  is 
proportional  to  ..’yr  .  This  explains  the  above  phenomenon 


that  when 


is  large  if  Mm  is  also  large  then  the  effect 


expansion  distance  is  short.  Comparing  to  the  limiting  number 
of  V  for  the  condition  that  only  turning  but  no  accelera- 

r . 

tion  occurs,  if  we  consider  that  the  same  magnitude  of  limiting 
number  should  be  applied  to  the  second  order  partial  derivative 
for  the  condition  that  only  acceleration  but  no  turning  occurs, 
then  approximately  we  get  < -  •;*  or 

\  •  s  (  T  2  ' 

V-  "  ' .  J  i  *  •  Otherwise  the  calculated  X  and  X  will  be  too 

•  s  p 

large.  Of  course,  if  it  turns  as  well  as  accelerates,  then 
both  limiting  regions  will  have  to  be  even  more  narrow 


However,  it  is  not  a  linear  additive  relationship.  Using  the 
Prandt 1-Meyer  flow  with  a  known  rigorous  solution  to  verify 
this  point,  we  found  that  this  rule  is  still  feasible.  The 
relation  between  and  M  is  shown  in  Figure  3  (the 

'  j .  >  t '  t 

figure  shows  the  condition  under  which  =  1 . 4  and  for  the 
generally  used  V  value  the  variation  of  this  curve  is  very 

small).  It  can  be  seen  that  for  the  same  .  ,  the  effec- 

v  u>  •• 

tive  expansion  distance  when  Mm  -  1.2  is  only  half  of  that 
for  the  incompressible  case  (under  the  pure  acceleration  and 
no  turning  condition).  Therefore,  for  supersonic  flow,  if 
the  acceleration  (deceleration)  of  the  gas  flow  is  severe 
then  the  use  of  only  up  to  the  second  order  derivatives  in  the 
mean  stream  line  method  is  not  acceptable.  But,  for  an  ordi¬ 
nary  good  cascade  (at  least  for  a  turbine  cascade)  the  accer- 
eration  (deceleration)  in  the  transonic  section  is  relatively 
slow  in  order  to  obtain  the  good  characteristics.  Therefore, 
it  is  not  hindering  the  use  of  the  mean  stream  line  method 
which  only  uses  up  to  the  second  order  derivatives  in  the 
calculation.  Figure  1  and  2  shown  before  are  such  examples. 

The  above  only  analyzed  the  effect  of  the  curvature  and 
accleration  on  the  mean  stream  line  on  the  partial  derivatives 
of  velocity.  The  effect  on  the  partial  derivatives  of  the  flow 
density  is  not  mentioned.  This  is  because  the  two  are  actually 
consistent  through  the  equation  of  the  one-dimensional  flow. 

The  latter  is  more  complicated  to  analyze.  Since  the  one- 
dimensional  isoentropic  flow  only  has  one  Independent  variable, 
therefore  It  is  sufficient  to  analyze  the  derivatives  of 
velocity  alone.  In  addition,  the  above  analysis  was  mainly 
carried  out  using  physical  diagrams.  Mathematically  it  is 
better  to  determine  its  error  or  convergence  after  the  deri¬ 
vation  of  higher  order  terms  of  the  series.  But  it  is  very 
complicated  and  may  not  be  necessary  for  an  engineering  calcula¬ 
tion  method  such  as  the  mean  stream  line  method. 


Figure  3.  The  /!'■■■'.' ■?  ^  ^  curve  {T-i.k) 

From  the  above  analysis,  we  can  explain  the  following 
facts : 

1.  The  mean  stream  line  method  is  more  suitable  for 

r  3 1  ro] 

impulsive  turbine  cascade  than  for  the  reversed  type  ’  . 

This  is  because  that  the  latter,  with  exception  of  the  possi- 

S  '  . 

bility  of  having  a  larger  •' ,  definitely  has  a 

'Stm  1  * 

which  is  much  larger  than  that  of  the  former.  Therefore,  it 
corresponds  to  the  fact  that  the  latter  has  one  more  factor 
affecting  the  inaccuracy  of  the  latter  than  that  of  the  former. 

2.  Sometimes  (such  as  the  reversed  type  cascades)  the  60% 

stream  line  is  used  as  the  stream  line  selected  for  expansion 

on  both  sides  which  is  more  effective  than  choosing  the  mean 

r  q  i 

stream  lineL  J.  This  is  because  the  former  sometimes  is  even 
closer  to  the  geometric  center  of  the  flow  boundaries  than 
the  latter  which  makes  the  expansion  distance  Ay  on  either 
side  smaller.  Therefore,  the  expansion  is  even  more  accurate. 

n  t 


3.  Tn  genera  1 ,  ♦  ».  •  ...  .  ...  j  of 

the  flow  boundary  to  :  •  ..  .  ;  .  planar 

cascade  is  not  accura*  ■.  •  ...  Is 

already  too  large.  It  g..;-  . . .  advantages 

of  the  mean  stream  1 1 ; .  ••  ■  f  * !  o  center 

stream  line  which  has  tr.  •  ■  :  •  -  •  ;  •  ,  :  :.  distance  and 

the  A’  on  that  line  is  ale  *  han  those  on 

either  side. 


Figure  4 .  Cascade  pattern  nr.  i  its  experimental  and  calculated 
results.  Key:  1.  relative  coordinate  in  the  arc  direction, 

2.  experimental  value,  >.  taed  arc,  4.  inner  arc,  5.  mean 
stream  line  method. 

As  an  example  purposely  leslened  to  test  the  method. 
Figure  4  shows  the  calculate-,!  results  of  a  supersonic  com¬ 
pressor  cascade  (actually  In  the  experiment  acceleration 


existed  within  the  flow  boundaries).  The  experimental  results 

are  obtained  from  Reference  [10].  Although  this  cascade  has 

been  operating  at  M  »  1.6  ^  2.1,  yet  the  calculated  results 

(especially  in  the  front  half  of  the  flow  boundary)  still  have 

a  certain  extent  of  reliability.  In  the  latter  half  of  the 

flow  boundary  M  number  has  reached  above  2.0.  Therefore, 

the  A  and  A  obtained  from  the  calculation  are  both  on  the  high 
®  P 

side.  From  this  we  can  see  that  the  use  of  the  mean  stream  line 
method  to  calculate  ordinary  transonic  cascades  (M  ir  1 . 3  1 . 4  ) 

should  be  feasible. 

Therefore,  the  presently  available  relatively  simple  mean 
stream  line  method  can  be  directly  extended  to  the  transonic 
planar  cascades.  Since  the  nature  and  the  equations  of  this 
method  are  not  changed,  all  its  advantages  and  disadvantages 

[21  [51 

still  remain  the  sameL  J *  .  They  will  not  be  repeated  here. 

As  for  the  design  experience  and  technique,  it  is  also  not 
difficult  to  extrapolate  from  the  experience  acquired  in  the 
subsonic  situation  as  long  as  the  characteristics  such  as  q(M) 
has  a  maximum  at  M  -  1  are  noticed.  For  example  in  the  selec¬ 
tion  of  the  position  of  the  mean  stream  line  in  the  forward 
problem.  In  the  pure  supersonic  region  it  shou1d  be  chosen 
closer  to  the  inner  arc  than  the  geometric  center  line  which 
is  opposite  to  that  in  the  subsonic  condition. 

In  order  to  further  increase  the  accuracy  of  this  method, 

it  is  possible  toadd  higher  order  derivatives  as  the  correction. 

But  the  calculation  equation  will  become  much  more  complicated. 

It  also  may  not  be  able  to  expand  very  far.  It  may  not  even 

be  better  than  the  situation  that  a  shorter  expansion  distance 

is  first  chosen  in  the  selection  of  the  stream  line  expansion 

to  obtain  a  new  stream  line  and  then  gradually  expand  towards 

the  outside  to  be  more  effective.  However,  such  a  method  does 

not  allow  the  use  of  the  advantages  of  the  mean  stream  line 

HI 

method  .  Another  possible  way  to  solve  the  high  M  flow 


problem  Is  to  use  the  mean  stream  line  method  to  solve  the 
transonic  region.  After  reaching  the  pure  supersonic  region, 
we  will  change  to  the  uncomplicated  supersonic  characteristic 
line  method  to  go  on  calculating.  The  total  calculation  would 
not  be  too  complicated  either. 

V.  Conclusions  and  Projections 

Through  qualitative  analysis  and  actual  calculation,  this 

paper  proved  that  the  existing  mean  stream  line  method 

(especially  the  analytical  solution  which  eliminates  the  error 

in  numerical  differentiation)  can  be  used  to  Calculate  flow 

in  an  ordinary  transonic  planar  cascade  and  still  maintains  the 

r?l  fqi 

unique  characteristics  of  this  methodL  LJ  .  For  example, 
results  with  sufficient  accuracy  can  be  obtained  using  a  slide 
rule  within  a  couple  of  days  (in  the  reverse  problem).  This 
paper  also  discussed  the  limitations  of  this  method. 

The  extension  of  the  method  introduced  in  this  paper  to 
the  flow  in  a  cascade  on  a  revolving  surface  should  also  be 
feasible  as  indicated  by  the  precedent  ^ 

In  addition,  as  for  whether  there  is  a  more  generalized 
meaning  in  the  reduction  of  a  multi-dimensional  flow  problem  to 
an  one-dimensional  one  followed  by  the  use  of  the  characteris¬ 
tic  that  the  one-dimensional  flow  equation  does  not  have  too 
much  difficulty  in  the  transonic  region  to  calculate  the  flow 
on  a  chosen  stream  line  then  to  expand  outward  and  finally  to 
obtain  the  transonic  flow  solution  is  yet  to  be  studied. 


•Fan  Yee-Chien  of  Nanking  Turbine  Electric  Machinery  Factory 
and  L.S.  Dzung  of  BBC  in  Switzerland  have  worked  in  this  area. 
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APPENDIX  A* 


In  Appendix  A,  we  briefly  introduced  the  use  of  flow 
function  expansion  to  obtain  flow  boundary  and  the  analytical 
solution  of  the  planar  flow  which  is  expanded  in  the  normal 
direction  of  the  stream  line  using  the  mean  stream  line  method. 
It  also  explained  that  the  analytical  solution  of  the  mean 
stream  line  method  can  also  be  used  for  compressor  planar  cas¬ 
cades  with  higher  density. 

I.  The  Use  of  Flow  Function  To 
Solve  for  the  Flow  Boundary 

For  every  z  on  the  boundary  of  the  flow  in  the  cascade, 
the  flow  function  has  the  relation  ^  d ^  .  If  we  consider 

y  as  a  function  of  =  ^  L  ^ then 

(A-1) 

and 


(A-2) 


The  partial  derivatives  of  y  with  respect  to  ip  in  equation 
(A-12)  can  be  obtained  from  equation  (A-l)  as: 

'  i _ i 

-V~  v,  (A-3) 
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(A-A) 


*This  appendix  is  mainly  a  combination  of  the  graduation  theses 
of  some  of  the  students  (Liu  Yen-hau,  Yuen  Chie-Quin,  Lee  Foo- 
ming,  Wu  Wen-ching,  Yu  Tso-chu,  etc.)  of  the  class  of  59  in  the 
department  of  modern  mechanics  at  China  Science  and  Technology 
University.  Since  it  has  never  been  published,  it  is  included 
here . 
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Since  the  original  mean  stream  line  method  only  uses  up 
to  the  second  order  derivatives  to  provide  sufficient 

accuracy,  therefore  the  flow  functj.uh  expansion  is  also  used 
up  to  the  term.  Substituting  the  above  three  equations 

into  (A-2),  we  then  can  obtain  the  equation  of  the  flow  boundary 
which  can  be  solved  directly: 


H 


(A -6) 


where  Ay^  is  the  expansion  distance  when  the  flow  is  assumed 
to  be  a  pure  one-dimensional  flow. 


(A-7 ) 


Using  this  method  to  find  the  solution  for  the  example 
in  Reference  [3]  and  comparing  the  result  with  that  obtained 
originally  using  the  third  order  algebraic  equation  method, 
we  realized  that  the  solutions  obtained  using  these  two  methods 
are  very  close. 


II.  The  Method  Used  to  Expand  in  the  Normal 
Direction  of  the  Selected  Stream  Line 


The  original  mean  stream  line  involves  expansion  in  the 
y  direction  which  is  convenient  in  the  calculation  of  planar 
cascades.  If  we  establish  a  new  coordinate  system  (l,n)  at 
each  expansion  point  (z,y)  on  the  mean  stream  line  where  1  is 
the  tangential  direction  of  the  stream  line  and  n  is  its  nor¬ 
mal  direction  with  the  direction  defined  as  pointing  from  the 
center  of  curvature  toward  the  stream  line  (See  Figure  A-l), 

V 
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then  it  is  possible  to  expand  in  the  new  coordinate  n  direction 
for  each  point.  At  this  time  the  forms  of  the  partial  differ¬ 
ential  equations  in  Reference  [3]  are  still  valid  with  respect 
to  the  new  coordinate  systems.  The  only  thing  is  that  all  the 
B_  values  should  be  0  looking  from  the  new  coordinates.  In 
addition,  the  superscript  '  should  be  considered  as  the  partial 
derivative  of  1  rather  than  2.  The  two  derivatives  should 
have  the  following  relationships: 


4  *'L  M  *•’  *•»  > 


(A-8) 


( A— 9 ) 


Therefore  of  all  the  .  .  ,  and  B(z)  are  originally  known, 
then  it  is  possible  to  calculate  all  the  and  ■■■;  ,  based 
on  the  above  t.v:o  equations.  Subsequently,  they  can  be  sub¬ 
stituted  into  the  partial  differential  equations  in  the  n  dir¬ 
ection  to  find  the  solution.  Actually  since  the  B  values 

m 

are  always  zero  in  the  new  coordinate  system,  it  is  easier  to 
obtain  the  partial  differential  equations  with  respect  t.o  the 
n  direction  than  getting  those  in  the  y  direction.  It  is  shown 
as  the  following  (where  Q  is  the  radius  of  curvature  of  the 
selected  stream  line). 

( A-10 ) 

lA-ii ) 

,A -U) 

( A- 13) 
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In  the  determination  of 


t  h*'- 


to  use  the  same  flow  function  ex;  x:.. 
also  should  be  similar  to  that  ext 


.  :  .  *  ion 


C  A — 1  h ) 


line 


Figure  A-l.  1-n  Coordi-  Figure  A-2.  The  boundary  With  the 
nate  system.  Return  Flow  Massing  by. 

The  equation  which  determines  the  velocity  distribution 
then  becomes: 
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*  ^  ‘ 
\  r. 
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(A-15) 


If  the  normal  direction  expansion  method  is  used  to  solve 
the  planar  cascade  problem,  although  the  calculation  is  less 
In  the  determination  of  partial  derivatives,  yet,  because  there 
is  another  transformation  process  to  turn  the  flow  boundary  from 
the  (1,  n)  coordinate  system  back  to  the  original  (Z,  y) 
coordinate  system,  the  total  amount  of  calculation  is  not  going 
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to  be  too  much  less.  Furthermore,  the  normal  direction 
expansion  method  makes  it  most  difficult  to  calculate  the 
velocity  distribution  in  the  most  important  inlet  and  outlet 
regions  on  the  back  arc.  Therefore,  it  may  not  be  suitable 
to  use  a  pure  normal  direction  expansion  method.  However, 
for  some  flow  boundaries  where  the  cross  sections  of  the 
inlet  and  outlet  are  perpendicular  to  the  mean  stream  line  or 
where  there  is  a  large  turn,  it  is  better  to  use  this  method. 

For  example,  the  boundary  where  the  return  flow  passes  by  in  a 
return  flow  turbine  (See  Figure  A-2)  is  such  a  case.  We  have 
used  this  method  to  design  the  flow  boundary  of  a  natural  gas 
expansion  turbine.  The  preliminary  result  was  satisfactory. 

III.  The  Feasibility  of  Using  the  Mean  Stream 
Line  Method  to  Calculate  Axial  Flow 
Compressor  Planar  Cascades 

Reference  [3]  pointed  out  that  it  is  possible  to  use  the 
existing  mean  stream  line  method  to  design  planar  cascades  of 
the  axial  flow  compressors.  In  order  to  verify  this  view, 
we  have  actually  solved  a  series  of  forward  and  reverse  prob¬ 
lems  of  this  type  of  cascades.  The  forward  problem  is  to 
find  solutions  for  cascades  with  known  velocity  distributions 
such  as  the  NACA  65  series  and  the  NGTE  C  series  cascades. 

The  range  of  variation  of  the  cascade  density  calculated  is 
1.06  ~1.20.  The  range  of  variation  in  the  geometric  deflection 
angle  of  the  blade  profile  is  24.6°  65°.  The  results  obtained 

from  the  forward  problem  are  more  consistent  to  those  obtained 
experimentally.  This  also  means  that  for  axial  flow  com¬ 
pressor  planar  cascades  with  higher  density,  it  is  still  possible 
to  find  the  solution  using  the  existing  mean  stream  line  method 
and  the  first  three  terms  of  the  Taylor  series.  It  breaks  the 
tradition  that  such  a  method  was  only  able  to  be  used  in  the 
calculation  of  turbine  cascades. 


APPENDIX  B 


Appendix  B  gives  all  the  expressions  and  numerical 
tables  of  the  5^(6)  functions  in  equations  (1)  -  (4)  to  be 
used  during  the  course. 
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THEORY  OF  THREE-DIMENSIONAL  FLOW  IN  TRANSONIC  TURBINE  MACHINERY 
WITH  SHOCK  WAVE  INTERRUPTIONS 

Hsu  Tsen-chung 

ABSTRACT 

In  order  to  be  consistent  with  the  three-dimensional  spatial 
flow  in  turbine  machinery  with  shock  waves,  this  paper  derives 
the  "four-dimensional"  basic  equations  of  the  S^  and  flow  sur¬ 
faces  through  the  extension  of  the  flow  surface  concept  to  the 
unsteady  flow  condition.  They  are  also  applicable  to  spatial 
three-dimensional  unsteady  flow.  Using  the  theory  of  the  charac¬ 
teristic  line  for  these  equations,  we  can  obtain  the  mutual  rela¬ 
tionships  between  the  characteristics.  From  them,  we  can  deter¬ 
mine  the  boundary  conditions  of  the  steady  flow.  These  basic  equa¬ 
tions,  boundary  conditions  and  applicable  initial  conditions  can 
completely  establish  the  transonic  flow  problem  with  shock  waves 
for  both  types  of  the  flow  surfaces.  Thus,  based  on  the  procedures 
recommended  in  this  paper  to  obtain  the  complete  solution  of  the 
three-dimensional  flow  and  by  choosing  the  proper  difference  equa¬ 
tions,  the  numerical  solution  can  then  be  obtained. 

I.  INTRODUCTION 

In  order  to  obtain  the  solution  of  three-dimensional  flow  in 
subsonic  and  supersonic  turbine  machinery,  [1]  presented  the  and 
S2  relative  flow  surface  theory  to  establish  the  basic  equations 
and  boundary  conditions  on  these  two  types  of  flow  surfaces  and 
also  presented  procedures  to  solve  these  equations.  For  over  20 
years,  this  theory  has  been  widely  used.  This  practice  indicates 
that  this  two  flow  surface  thoery  is  not  only  the  calculation  method 
for  three-dimensional  flow  in  the  turbine  machinery,  but  also  the 
theoretical  basis  of  the  design  method  of  such  machinery.  The  cal¬ 
culation  and  iteration  processes  of  the  two  types  of  flow  surfaces 

tf 


are  not  only  the  procedures  to  solve  the  three-dimensional  flow 
problem,  but  also  those  to  design  three-dimensional  flow  turbine 
machinery . 

Recently,  with  the  development  of  transonic  turbine  machin¬ 
ery,  there  has  been  a  desperate  need  to  calculate  three-dimensional 
flow  with  shock  wave  interruptions.  In  [2],  it  was  pointed  out, 
after  the  establishment  of  the  basic  equations  of  three-dimen¬ 
sional  flow  and  the  integral  form  on  the  two  kinds  of  flow  surfaces 
and  the  derivation  of  the  relations  of  parameters  in  front  and 
behind  the  shock  wave  under  various  conditions,  that  in  order  to 
solve  the  three-dimensional  flow  problem  we  must  calculate  a  two- 
dimensional  flow  which  is  consistent  with  the  three-dimensional 
flow  on  each  of  the  two  surfaces,  and  carry  out  iterations  between 
the  calculations  of  the  two  surfaces. 

The  establishment  of  this  type  of  two-dimensional  flow  is  a 
fixed  solution  problem.  Based  on  the  results  in  [2],  it  cannot  be 
carried  out  on  an  assigned  flow  surface.  It  must  extend  the  steady 
flow  surface  model  presented  in  [1]  to  the  unsteady  flow  condition* 
It  transforms  three-dimensional  flow  to  the  two  types  of  flow  sur- 
surfaces  which  vary  with  time.  This  paper  first  briefly  introduces 
this  unsteady  surface  and  the  geometric  relations  on  it.  It  also 
shows  the  relation  between  three-dimensional  flow  surfaces  at 
every  instance. 

It  is  then  followed.  In  the  second  an  1  third  sections  of  this 
paper,  by  the  establishment  of  the  basic  equations  on  the  two  flow 
surfaces  which  are  consistent  with  the  three-dimensional  flow.  They 
are  applicable  in  the  solution  of  the  unsteady  flow  and  the  cal¬ 
culation  of  steady  flow  with  shock  waves. 

1 - 

Comrade  Huang  Ray  Tsien  has  pointed  out  that  it  is  possible  to 
extend  the  steady  flow  surface  to  the  unsteady  flow  surface. 
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In  the  fourth  section  of  this  paper,  the  characteristic  rela¬ 
tions  on  the  two  types  of  surfaces  are  derived  on  the  basis  of 
characteristic  line  theory.  It  presents  the  boundary  conditions 
at  upstream  and  downstream  boundaries  and  on  the  flow  surfaces 
which  are  suitable  for  solving  the  problem  of  steady  flow  with  shock 
waves.  The  initial  conditions  can  be  given  using  the  usual  method. 
Thus,  this  specific  solution  problem  can  be  completely  established. 

Finally,  in  the  fifth  section  of  this  paper  we  briefly  des¬ 
cribe  the  procedures  used  to  solve  the  three-dimensional  flow  prob¬ 
lem  by  iteration  of  calculations  between  the  two  types  of  unsteady 
flow  surfaces. 

II.  UNSTEADY  FLOW  SURFACE  AND  ITS  GEOMETRIC  RELATION 

The  combination  of  all  the  stream  lines  of  all  the  points  on 
a  curve  in  the  flow  field  which  is  a  stream  line  forms  the  spatial 
flow  surface. 

In  unsteady  flow,  the  flow  surface  varies  with  time.  Therefore 
In  the  four-dimensional  space  including  time,  the  flow  surface  can 
be  expressed  as:  _ 

if.lfc  ■</>  &ij.=pr. 

Thus 


The  unit  vector  in  the  normal  direction  N(.Nr,Nu  ,Nz  ,Nt ) 
satisfies 


Therefore,  on  the  four-dimensional  flow  surface: 


(1) 


'tfdr+A' niyi-Aj d y  +A(.  U<d ir=*o 

:  f/t  _  i  _/.  '  J ~t  ;  :A  "t 
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On  the  other  hand,  there  is  a  flow  surface  at  any  instance 
in  the  space 

Its  unit  vector  n(nr,  n^,  )  is  a  function  of  time  t  and  it 

satisfies 


In  the  meantime,  from  the  definition  of  the  flow  surface,  for 
every  instance  we  have  _  _ _ 

71 -W  =t0. 

which  is 

"t*  \\X  Q 

It  should  be  pointed  out  that  the  corresponding  components  of 

the  three-dimensional  and  four-dimensional  unit  normal  direction 
-►  — — 

vectors  n  and  N  are  not  equal,  but  there  are  the  following  rela¬ 
tions  : 

Me 

~ ?lr  *~  — /V£ 7lr-J Au. - - ; —  (2) 

-  '  ,  .  '  ” 

With  the  above  relations  on  the  unsteady  flow  surface,  we  can 
transform  the  basic  equations  of  the  non-viscous  gas  in  a  relative 
coordinate  system  which  is  rotating  at  a  constant  angular  velocity 
as  shown  in  [1]: 


separately  to  the  two  types  of  unsteady  flow  surface.  It  should 
be  pointed  out  that,  as  a  closed  series  of  equations,  the  equa¬ 
tion  of  the  state  of  the  gas  and  other  relations  should  be 
Included  in  addition  to  the  above  equations. 

C8) 

and 

Because  these  algebraic  equations  maintain  their  forms  on 
the  flow  surface,  they  will  not  be  described  further  in  the  follow¬ 
ing  sections. 

III.  BASIC  EQUATIONS  ON  THE  UNSTEADY  FLOW  SURFACE  5 

Let  us  consider  the  unsteady  flow  surface  formed  by  the  stream 
lines  passing  through  a  certain  circle  located  upstream  of  the 
cascade . 
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Equations  (13)— (15)  can  also  be  expressed  using  I  and  S: 
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If  we  are  considering  a  steady  flow,  the_ f low  surface  does 
not  vary  with  time.  In  the  above  equations  JJ-== 1  £>  <p 

Then  these  equations  can  be  transformed  into  the  basic"  *  equations 
of  the  steady  flow  surface  corresponding  to  Equations  (3#b),  (Jio), 


c^l) ,  (39a),  (39b)  and  (39c)  in  [1].  In  other  words,  comparing 

to  the  equations  on  the  steady  flow  surface,  in  the  eouations 

on  the  unsteady  flow  surface,  besides  the  additional  fLil-l 
,,//  x?  '  +  .  hr 

and  —r—^ir  in  the  expressions  of  c  and  f  respectively.  In  the 
cf 

energy  equation,  there  is  a  similar  J  term. 


We  found  that,  different  from  the  basic  equation  of  transonic 
flow  with  shock  wave  interruptions  on  the  flow  surface  used  to 
date,  in  Equations  (12)-(l6)  b  varies  with  time  and  appears  in  the 

~~  term.  Simultaneously,  it  shows  the  effect  of  unsteady 
.  ?T  _ 

factors  in  f  and  J.  These  are  because  of  the  use  of  the  four¬ 
dimensional  steady  flow  surface  model  in  the  derivation  of  these 
equations  which  are  consistent  with  the  spatial  shock  wave. 

Considering  the  process  of  finding  for  a  complete  numerical 
solution,  it  is  more  accurate  to  carry  out  the  calculation  of  the 
difference  format  in  the  invariance  form  during  the  interruption 
due  to  the  shock  waves.  We  can  transform  Equations  (12)~(l6)  Into 
the  following  divergent  forms  (or  still  called  the  Invariance 
forms ) 

j&f)|  I  aVW,  W)  |  bTn,  f  ?U  b£ssso  _ !  _ _ !• _ I _ i_  _j  ;  )Qt) 
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or 


where 


IV.  BASIC  EQUATION  ON  THE  UNSTEADY  FLOW  SURFACE  S2 

Let  us  consider  the  unsteady  flow  surface  formed  by  stream 
lines  passing  through  a  radial  line  located  upstream  of  the  cascade 

•  5.(0  3,  It).  o__ 

Any  arbitrary  function  q  on  it  can  be  expressed  as 

f,Jt)  : 

By  taking  Equation  (2)  into  consideration,  the  partial  deri¬ 
vatives  of  q  on  the  four-dimensional  flow  surface  can  be  trans¬ 
formed  as: 
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where 


From  these  relations.  Equation  (3)  can  be  expressed  as 
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by  introducing 
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the  continuity  equation  can  be  rewritten  as: 
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Similarly,  Equations  (^)-(7)  are  as  the  following  respectively 
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If  we  use  I  and  S  to  express  it,  then  the  momentum  equation 
can  be  transformed  into: 
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Similarly  to  the  situation  on  the  S1  flow  surface,  in  the 
steady  flow  the  above  equations  can  be  transformed  into  equations 
(100),  (98),  C99a),  (96a),  (96b)  and  (96c)  in  [1],  respectively. 
The  difference  between  the  steady  and  unsteady  S2  flow  surfaces  is 
also  the  same  as  the  situation  for  the  S-^  flow  surface. 


Equations  C30)— (3^ )  can  be  transformed  into  the  divergent  forms 


or 


where 
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It  should  be  pointed  out  that  the  basic  equations  on  the  two 
types  of  unsteady  flow  surfaces  are  applicable  in  slowing  the 
unsteady  flow  as  well  as  the  steady  flow  with  shock  waves.  The 
difference  in  the  two  specific  solution  problems  is  in  their 
boundary  conditions. 

V.  BOUNDARY  CONDITIONS  AND  INITIAL  CONDITIONS 

In  the  two  previous  sections,  we  have  already  established  the 
hyperbolic  partial  differential  equations  on  the  two  types  of 
unsteady  flow  surfaces.  The  following  is  a  presentation  of  the 
corresponding  boundary  conditions  and  initial  conditions. 

In  the  calculation  of  the  initial  value-boundary  value  problems, 
the  parameters  on  the  surfaces  and  at  the  upstream  and  downstream 
boundaries  are  related  to  the  parameters  of  the  flow  a  moment  before. 
Thus,  we  should  use  the  characteristic  line  theory  [3]  of  the  hyper¬ 
bolic  equations  to  derive  the  relative  relations  between  the  char¬ 
acteristics  in  order  to  determine  the  boundary  conditions  [**,5,6,73, 

Cl )  The  boundary  condition  on  the  S.  flow  surface 

At  this  time  we  can  obtain  the  following  relations: 
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where  y  represents  the  angle  between  the  outer  normal  line  and  the 
z  axis  at  a  certain  point  on  the  boundary.  Equations  (42)  and  (43) 
are  the  flow  characteristic  relations  while  Equations  (44)  and  (45) 
correspond  to  the  wave  characteristic  relations  of  gas  propagating 
at  a  speed  -a  and  a,  respectively. 

For  the  upstream  boundary:  Let  us  assume  that  ^,0  ^  ^  ;  then 

only  Equation  (44)  can  be  used.  Therefore,  three  parameters  must 
be  given.  The  other  parameters  can  be  obtained  from  Equation  (44). 
Please  notice  that  because  y  =  ir  at  this  time,  and  the  three  given 
parameters  (e.g.,  pQ,  pQ  and  gas  flow  angle)  do  not  vary  with  time 
and  the  parameters  along  the  ^  direction  are .homogeneous ,  Equa¬ 
tion  (44)  can  be  drastically  simplified. 


For  the  downstream  boundary:  If  the  exit  gas  flow  is  subsonic 
at  this  time,  then  Equations  (42),  (43)  and  (45)  are  usable  rela¬ 
tions.  It  still  needs  one  more  given  parameter  (such  as  the  re¬ 
verse  pressure).  Please  note  that  y  =  0  at  this  time  and,  the 
above  three  equations  can  also  be  simplified.  If  the  exit  gas  flow 
is  supersonic,  then  all  four  relations  can  be  used.  We  can  no 
longer  provide  any  given  parameter. 


For  the  blade  surfnoo:  Equations  (42),  (43)  and  (45)  are 
all  applicable.  The  required  boundary  condition  is  that  the  flow 
velocity  is  tangential  to  the  surface. 

In  addition,  there  is  a  periodic  condition  in  the  calculation 
for  the  flow  surface.  Its  calculation  can  be  carried  out  based 
on  the  same  method  as  that  for  the  internal  points.  Cf  course, 
when  the  shock  wave  extends  to  the  outside  of  the  duct,  whether 
to  choose  the  periodic  condition  of  the  suction  surface  in  the 
pressure  surface  requires  a  specific  analysis. 

(II)  Boundary  conditions  on  the  S~  flow  surface 


The  following  relations  can  be  obtained  at  this  time 

Tn<  ,  ,J*{  //3v  Tivi.  ,Ii 4  ,  /IPf 

s7(jr+"'' rr  >3  *f$J~ 
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where  y*  represents  the  angle  between  the  outer  normal  direction 
of  a  point  on  the  boundary  and  the  z-axis.  Equations  06)  and 
(47)  are  the  flow  characteristic  equations,  while  Equations  (48) 
and  (49)  are  the  wave  characteristic  relations  of  the  correspond¬ 
ing  gas  propagating  at  a  speed  of  -a  and  a,  respectively. 

For  the  upstream  boundary:  Let  us  assume  0'  0.  >  then 

only  Equation  (48)  can  be  used.  Therefore,  it  needs  three  given 
parameters.  The  other  parameter  can  be  obtained  from  Equation  (48). 
Please  note  that  because  y'  =  tt  and  the  three  given  parameters 
(e.g.,  Pq,  and  M  number)  are  time  independent  and  all  the  para¬ 
meters  in  the  r  direction  are  homogeneous.  Equation  (48)  can  be 
greatly  simplified. 

For  the  downstream  boundary :  If  the  exit  flow  is  subsonic  at 
this  time,  then  Equations  (46),  (47)  and  (49)  are  the  applicable 
relations.  It  still  requires  one  given  parameter  (such  as  the 
reverse  pressure  or  Vu  r).  Please  note  that  at  this  time  y'  =  0, 
the  above  three  equations  can  also  be  greatly  simplified.  If  the 
exit  gas  flow  Is  supersonic,  then  all  four  relations  can  all  be 
used.  It  cannot  be  provided  with  any  given  parameters. 

For  the  up  and  down  wails:  Equations  (46),  (47)  and  (49)  are 
all  applicable.  The  given  boundary  condition  is  that  the  flow 
velocity  is  in  the  tangential  direction  to  the  wall  surface. 


(III).  Initial  conditions 

It  has  already  been  proven  [8]  by  the  theory  that  the  choice 
of  the  initial  conditions  does  not  affect  the  final  steady  solu¬ 
tion.  However,  if  the  selection  of  the  initial  conditions  is  very 
improper,  it  may  cause  excessive  oscillations  and  damge  the  stab¬ 
ility  of  the  numerical  solution.  Therefore,  the  initial  values 
should  be  chosen  so  that  there  is  no  sudden  variation  for  all  the 


flow  parameters .  In  the  meantime,  the  selection  of  the  initial 
values  also  affect  the  speed  reaching  a  steady  solution.  Con¬ 
sidering  these  factors,  we  can  choose  the  solution  of  the  sub¬ 
sonic  flow  or  the  parameter  distribution  which  satisfies  the 
continuity  equation  on  the  geometric  center  line  and  appears  to 
be  smooth  elsewhere  as  the  initial  condition. 

VI.  THE  PROCEDURES  OF  OBTAINING  THE  COMPLETE  SOLUTION  OF  THREE- 

DIMENSIONAL  FLOW 

We  have  already  transformed  the  problem  of  the  three-dimen¬ 
sional  flow  with  shock  waves  in  turbine  machinery  into  two  families 
of  initial  value-boundary  value  problems  of  hyperbolic  equations 
with  three  variables.  It  also  presents  the  methods  to  determine 
the  initial  and  boundary  conditions  and  then  completely  establishes 
the  specific  solution  problems  for  the  two  types  of  unsteady  four¬ 
dimensional  flow  surfaces.  In  order  to  obtain  the  solution  of  the 
three-dimensional  flow,  besides  the  steady  solutions  of  the  initial 
value-boundary  values  problems  on  the  two  families  of  flow  surfaces, 
we  must  carry  out  iteration  between  each  other  until  it  becomes 
convergent . 

In  order  to  save  the  calculation  time,  we  can  use  the  time  pro¬ 
gressing  and  flow  surface  iteration  method.  For  the  calculation  of 
each  step  of  Atj  or  several  steps  of  EAt  j  ,  we  carry  out  an  iteration 
between  the  two  families  of  flow  surfaces.  This  involves  giving 
the  results  of  each  flow  surface  calculation  and  the  shape  of  each 
four-dimensional  flow  surface  (see  Equation  (1))  to  the  other  flow 
surface  to  carry  out  the  calculation  for  the  next  time  interval. 

It  repeats  itself  regardless  of  whether  a  steady  solution  is 
reached  in  each  calculation.  When  the  steady  solution  is  nearly 
obtained,  iteration  between  flow  surfaces  is  required  for  every 
time  interval.  Of  course,  it  is  not  required  to  have  a  completely 
steady  solution.  Finally,  when  steady  solutions  are  obtained  on 
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all  the  flow  surfaces,  the  final  solution  of  the  entire  three- 
dimensional  flow  is  obtained. 

Such  procedures  to  obtain  the  solution  clearly  demonstrate 
what  has  been  pointed  out  in  [2]  that  this  method  cannot  be  used 
to  solve  a  transonic  flow  with  shock  wave  on  a  given  flow  surface. 
In  other  words ,  it  Is  through  the  continuous  adjustment  of  the 
position  and  shape  of  the  flow  surface  during  the  solution  seeking 
process  which  makes  it  possible  to  obtain  the  two  dimensional  shock 
waves  which  are  consistent  with  the  spatial  shock  wave  on  the  two 
families  of  flow  surfaces.  The  above  procedures  closely  combine 
the  calculated  shock  wave  on  a  single  flow  surface  with  the  spatial 
shock  wave  so  that  the  time  needed  to  obtain  a  steady  solution  can 
be  reduced. 

VII  CONCLUSIONS 

On  the  basis  of  extending  the  three-dimensional  flow  surface 
model  to  the  four-dimensional  unsteady  flow  surfaces,  this  paper 
transformed  the  ordinary  equations  of  the  non-viscous  adiabatic  gas 
flow  in  turbine  machinery  to  two  types  of  unsteady  flow  surfaces 
S-^  and  S2  and  then  re-established  the  corresponding  basic  equations 
Comparing  these  equations  with  those  for  corresponding  steady  flow 
surfaces,  we  found  that,  in  addition  to  the  presence  of  the  partial 
derivatives  of  the  thickness  of  the  flow  sheet  vs.  time,  the  shape 
of  the  four-dimensional  flow  surface  has  a  certain  contribution 
toward  c  (or  c')  and  the  "flow  sheet  force".  In  the  energy  equa¬ 
tion,  there  is  also  a  related  term  e(or  e').  These  equations  are 
applicable  to  the  unsteady  flow  or  the  steady  flow  with  shock  waves 

It  should  be  pointed  out  that  the  equations  used  to  obtain 
solutions  to  the  transonic  flow  problem  in  the  literature  to  date 
on  the  flow  surfaces  (mainly  revolving  surfaces)  are  for  time  inde¬ 
pendent  flow  surfaces.  The  thickness  of  the  flow  sheet  is  also  not 


related  to  time.  Therefore,  the  shock  wave  thus  calculated  is 
not  a  part  of  the  real  shock  wave  in  the  space.  Actually,  the 
basic  equations  used  in  the  literature  are  not  obtained  from  the 
basic  equations  used  to  solve  for  the  spatial  shock  waves.  They 
are  merely  equations  with  the  arbitrary  addition  of  unsteady  terms 
to  the  steady  flow  equations.  The  result  is  that  it  failed  to  con¬ 
sider  the  close  relations  between  the  two  families  of  flow  sur¬ 
faces  under  the  conditions  that  shock  wave  existed. 

Based  on  the  characteristic  curve  theory,  we  derived  the  char¬ 
acteristic  relations  between  the  two  types  of  surfaces  and  further 
determined  the  boundary  conditions  for  the  condition  of  steady 
flow  with  shock  waves.  There  is  a  certain  degree  of  flexibility  in 
the  selection  of  the  initial  conditions.  The  proper  choice  of  the 
initial  condition  can  completely  define  the  specific  solution  prob¬ 
lem  of  transonic  flow  with  shock  wave  on  the  two  types  of  flow 
surfaces . 


Considering  the  close  relation  between  the  flow  surfaces  when 
there  are  shock  waves,  in  order  to  obtain  the  solution  of  the  three 
dimensional  flow,  we  can  use  the  time  progressing  and  flow  surface 
iterations  method  to  save  the  calculation  time.  Thus,  after  the 
difference  format  is  chosen,  we  can  obtain  the  solution  numerically 
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